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ABSTRACT 

Close-in giant planets are thought to have formed in the cold outer regions of planetary systems and 
migrated inward, passing through the orbital parameter space occupied by the terrestrial planets in 
our own Solar System. We present dynamical simulations of the effects of a migrating giant planet on a 
disk of protoplanetary material and the subsequent evolution of the planetary system. We numerically 
investigate the dynamics of post-migration planetary systems over 200 million years using models with 
a single migrating giant planet, one migrating and one non-migrating giant planet, and excluding the 
effects of a gas disk. Material that is shepherded in front of the migrating giant planet by moving 
mean motion resonances accretes into " hot Earths" , but survival of these bodies is strongly dependent 
on dynamical damping. Furthermore, a significant amount of material scattered outward by the giant 
planet survives in highly excited orbits; the orbits of these scattered bodies are then damped by gas 
drag and dynamical friction over the remaining accretion time. In all simulations Earth-mass planets 
accrete on approximately 100 Myr timescales, often with orbits in the Habitable Zone. These planets 
range in mass and water content, with both quantities increasing with the presence of a gas disk 
and decreasing with the presence of an outer giant planet. We use scaling arguments and previous 
results to derive a simple recipe that constrains which giant planet systems are able to form and 
harbor Earth-like planets in the Habitable Zone, demonstrating that roughly one third of the known 
planetary systems are potentially habitable. 

Subject headings: astrobiology - methods: N-body simulations — planetary systems: formation 



1. INTRODUCTION 

More than 2 00 giant planets are kno wn to orbit main- 
sequence stars dButler et al.l (|2006T ) : see lSchneiderl (2006) 
for recent results) and all but 5 have semi-major axes 
within the orbit of Jupiter; more than half of known plan- 
ets reside within 1 AU of their parent star. There are also 
a surprising number of planets at very small semi-major 
axes: 22% of currently known extrasolar planets have or- 
bital radii less than 0.1 AU, and 16% are located within 
0.05 AU of the central star. Limitations on the current 
observational techniques do not allow a complete sample 
of planets aroun d solar-type stars beyond appro ximately 
3 AU (see, e.g., iTabachnik fc Tremaind ([20021 ) L but it 
is clear that there is a significant population of plane- 
tary systems with giant planets at small orbital radii. 
The occurrence of close-in giant planets is surprising be- 
cause theoretical models predict that the formation of 
giant planets in the hot inner regions of a protoplan- 
etary disk would be dif ficult under realistic conditions 
(|Bodenheimer et al.l l2000h . To explain the observations 
of massive planets at very small distances from their 
parent stars, a mechanism is necessary to move them 
from where they formed to where they currently reside; 
this process is c ommonly known as planetary migration 
(|Lin et al.lll996h . With the inclusion of migration of a 
giant planet to small orbital radii, theories on the evo- 
lution of solid bodies in the inner disk developed for our 
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own planetary system must be re-examined. 

1.1. Giant Planet Formation and Migration 

Favored theories of giant planet formation center 
around two main paradigms, commonly called the core 
accretion model and the gravitation al instability model . 
The bottom-up core accretion model (|Pollack et ailll996l : 
lAlibert et alJl2005t iHubickvi et al.ll2005fi " requires the ac- 
cretion of planctesimals into a solid core of ~ 5-10 M®, 
massive enough to initiate runaway gravitational infall 
of a large gaseous envelope. The formation of this mas- 
sive core requires a high surface density of solids, which 
is difficult to achieve in the hot inner disk. A jump 
in the surface density is believed to occur just past the 
"snow line" , where the disk temperature drops below the 
freezin g point of wat er (at ~ 170 K in protoplanetary 
disks ( Hayashi! 119811 ) ) and ice is available as a build- 
ing block. Isolation masses also increase with orbital 
distance r for surface density profiles flatter than r~ 2 
(|Lissauerl H987). also favoring formation at larger orbital 
distances. The top-down gravitational instability mode l 
(|Bosslll998l 120001: iMaver et al1l2002t iDurisen et ai1l2005l) 
assumes a gaseous protoplanetary disk massive and cold 
enough to become gravitationally unstable to collapse, 
resulting in a massive gaseous planet with little or no 
solid core. For realistic disk masses, this process only 
becomes efficient a t orbital radii greater than ~10 AU 
(jMaver et all 120021) . Both are believed to be viable un- 
der certain physical conditions, but it is still unclear how 
each mechanism (or combination of mechanisms) func- 
tions in realistic circumstellar environments; even cal- 
culations of the core masses for the gia nt planets in 
our own system have la rge uncertainties (fGu illot 1999; 
iSaumon fc Gui llot 2004). Regardless of which mecha- 
nism forms planets in extrasolar planetary systems, it is 
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clear that giant planets detected as small orbital radii 
were unlikely to form at their present locations. 

A variety of potential migration mechanisms 
have been proposed to operate on massi ve plan- 
ets, including planet-planet s catter i ng (IFarinella 
19801: IWeidenschilling fe Marzari) 119961: iRasio fc Ford 



1996 ) , planetesima l scattering ( Fernandez fe Id Il984t 



Murray et alj 119981). and a several types of gas-p lanet 
interactions (see iPapaloizou fc Terqueml ([20061 ) for 



a review 



The most simple and robust hypothesis 
is known as Type II m i gration. Type II migration 
(|Papaloizou fc Lid 11984 iLin et all 1 19961 ) works on 
planetary bodies that are massive enough to open an 
azimuthal and vertical gap in the gas disk, locking 
the planet into common orbital evolution with the gas 
disk. Viscosity in the gaseous disk, thought to be due 
primarily to th e magneto-rotational instability (MRI, 
iGammiel (|1996l) ). results in orbital decay and infall of 
the gas onto the central star (as e videnced by accretion 



and stellar a ctivity in young stars (iMuzerolle et al.l L 
ICalvet et ail 12004 lEisner et alj 120051 )). Thus, the giant 
planet loses angular momentum and migrates inwards 
coupled to the gas disk. Simulations suggest Type II 
migration timescales range between 10 5 to a few times 
10 5 years, depending o n the di sk and planet co n dition s 
([Nelson fc Papaloizoul 12004 iD'Angelo efafl 120031 ) . 
Planets halted at very small radii most likely ceased 
migration due to gas evacuation around the central star 
or planetary mass loss onto the c entral star, creating the 
"hot Jupiter" population (see IPapaloizou fc Terqueml 
(2006) for further explanation). However, for orbital 
radii beyond approximately 0.1 AU these stopping 
mechanisms would not be effective, and other expla- 
nations are required for extrasolar planets detected in 
orbits between 0.1 and a few AU. The intermediate 
stopping distances may be a result of dissipation of 
the gas disk during migration ([Trilling et al.lTl99l or 
a "d ead zone" where the MRI is ineffective ( Gammiel 
1996), but these mechanisms rely on extreme fine-tuning 
of parameters or questionable physical conditions to 
explain the distribution of orbital radii for extrasolar 
planets. 

1.2. Terrestrial Planet Formation 

Standard theories of the evolution of a planetary 
system suggest that a circumstellar disk will proceed 
through various accretionary stages culminating in the 
final architecture of a stable planetary configuration. Co- 
agulation of the inceptive dust particles occurs through 
collisional sticking until the solid bodies become mas- 
sive enough to decouple from the surrounding gas disk 
and settle to the disk midplane, resul ting in meter- 
sized objects on time scales of 10 4 years ([Lissauerl 119931 : 
iBeckwith et al.1 120001) . In this si ze range orbital decay 
due to gas drag is very rapid (jWeidenschillingJ 119771 ) 
and accretion must happen very quickly to avoid in- 
fall onto the central star and overly extensive mixing. 
The accretion rate can be enhanced by differential mi- 
gration r ates ( Wcidcnschill ingl [19771). gravitational in - 
stability (|Goldreich fc Wardlll973l : lYoudin fc Shdl2002l) . 
and/or concentration of meter-sized bodies in spiral den- 
sity waves generated in the gaseous disk ([Rice et al.l 
2004). Once the largest bodies reach ~1 km in size, their 
gravitational cross-section becomes larger than their ge- 



ometric cross-section and they begin the phase known as 
"runaway growth". During this phase the largest bod- 
ies grow faster than the smaller bodies due to gravita- 
tional focusing and damping due to dy namical friction, 
resul ti ng in a widening mas s dispe rsion (Grecnbcr g et all 
1978] IWetherill fc Stewart! 119891 : IWeidenschilling et all 
19971 ). "Runaway growth" transitions into "oligarchic 



growth" when the velocity dispersion of planetesimals 
becomes comparable to the e scape speed of the largest 
embryos (|Kokubo fc Idalll998l ), stalling accretion for the 
largest bodies and decreasing the embryo mass disper- 
sion. Oligarchic growth ends after approximately 10 6 ~ 7 
years when the oligarchs have depleted their "feed- 
ing zones" sufficiently such that dynamical friction is 
no lon ger effective, and embryos begin to scatter each 
other ([Kokubo fc Idal I1998L |2000[ ). The final "chaotic 
phase" of planet growth proceeds through scattering 
and collisions between the large protoplanets and fi- 
nal clearing of the remaining planetesimals to produce 
a stable planetary system after more than 10 8 years 
(IWetherilllll996l:IChamberdl20011 [Raymond et aLll2006d: 
iKenvon fc Bromlevll2006i r 

Though the evolutionary state of the protoplanetary 
disk at the inception of giant planet migration is un- 
certain, there are several reasons to believe that rel- 
atively large objects would exist in the inner disk by 
the time migration begins. Recent giant planet for- 
mation models assuming core accretion as the domi- 
nant mechanism give full forma tion timescales ranging 
from 5 Myr to less than 1 Mvr (iRice fc Armitagel 120031 : 
lAlibert et al.|[2005t iHubickvi et al.ll2005l ). corresponding 
to the oligarchic growth stage in standard terrestrial 
planet models. Observations of gas and dust in pri- 
mordial disks give d isk lifetimes of ~ 10 7 years for gas 
(|Bricefio et a l.1 [20011: lHaisch et alJl200lD and less than 3 
Myr for dust ( Sil verstone et all 120061 ) . suggesting that 
solid bodies form long before the gas disk dissipates. For 
our own Solar System, radioactive dating of meteorites 
give equivalent ages for both the earliest chondrule for- 
mation and differentiation in asteroidal bodies in our own 
Solar S ystem, supporting rapid evolution of large solid 
bodies (|Kleine et al.l 120021 ). Timescales for giant planet 
formation by gas instability may be as short as 10 3 years, 
but in the massive disks required for gas instability to 
function the migration mechanis ms may be halted until 
much of the gas has dissipated ([Boss! [20051 ) . Likewise, 
if massive disks are indeed required for giant planets to 
form quickly enough for migration to occur, formation 
timescales in other parts of the disk as would most likely 
be shortened as well. 

One of the most interesting problems in the forma- 
tion of planetary systems is the distribution of volatiles 
compared with refractory materials in different plane- 
tary bodies. Naturally, the final composition of planets 
formed through the scenario described above depends on 
both the initial distribution of material in the protoplan- 
etary nebula and the subsequent dynamical evolution of 
the system. The initial composition of solids in a cir- 
cumstellar disk is expected to f ollow a basic condensa- 
tion sequence (|Grossmad IT972). but the effects of ra- 
dial transport of grains and mm to meter-sized objects 
through turbulence or radial drift have the potential to 
dramatically re- arrange this early compositional gradient 
([Stevenson fc LunindTl98§ : IStepinski fc Valageasl fl997l : 
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ICiesla fc Cuzzil [20061 ). However, evidence from our own 
system suggests that late-stage protoplanetary material 
may maintain the gross characteristics of the basic con- 
densation sequence: analysis of chondrites in meteorites 
from parent bodies in different regions of the inner solar 
system suggest a constant increase from alm ost no wa- 
ter a t 1 AU to almost 10% water at 2.5 AU (|Abe et al.l 
2000), with cometary bodies thought to contain at least 
half of their mass in water ice . The dry nature of as- 
teroidal material in the vicinity of Earth further sug- 
gests that continued delivery of water-rich material must 
have been necessary to produce the c urrent water inven- 
tory, either through comet impacts (|Owen fc Bar-Nunl 
|1995|) or through accretion o f disrupted asteroidal ma- 
terial ( Morbidelli et al. 2000). The importance of radial 
transport of material in defining the final composition of 
terrestrial planets has been confirmed by N-body simu- 
lations, which suggest that late-stage accretion of water- 
rich material is largely influenced by the characteristics 
of giant plan ets in these systems, notably their orbital 
eccentricities ([Chambers fc Cass en 2002; Raymond et al.l 
[200llRavmondll2006D . 

1.3. Previous Work 

With the uncertainty inherent in the cornucopia of 
forces that may be functioning on objects in various dif- 
ferent size regimes in circumstellar disks, untangling the 
complex dynamics during the formation and evolution of 
solid bodies in the inner regions of a protoplanetary disk 
becomes difficult. Several studies concerning the hab- 
itability of planetary systems in a variety of conditions 
made the simple assumption that sys tems with close-in 
giant planets could n ot be habitable (|Lineweaverl 120011 : 
iGonzalez et al.l 120011 ). Give n the 'hot Jupite r '-stell ar 
metallicity corr elation (e.g., iFischer fc Valentil ([20051) ). 
this placed limits on the galactic locales likely to har- 
bor Ea rth-like planets (the so-c alled 'Galactic Habitable 
Zone'; iLineweaver et al.l (|200 4f)). More recent studies 
have concentrated on exploring the dynamical effects of 
giant plane t migration on a simplified disk, with vary- 
ing results. lArmitagg (|2003h . assuming planetesimals in 
the inner disk would be destroyed by a migrating giant 
planet, analyzed the evolution of the post-migration disk 
surface density with a simple dust coagulation formu- 
lation. They found that if migration begins after sig- 
nificant embryo formation has occurred, post-migration 
formation would be unlikely due t o dust depletion and in- 
suffici ent remaining material. In iMandell fc Sigurdssonl 
( 2003) we examined the assumption of planetesimal de- 
struction during migration by analyzing the probability 
of terrestrial-mass objects surviving the inward migra- 
tion of a giant planet. We used dynamical modeling to in- 
vestigate the survival rate and the eccentricity and semi- 
major axis distribution of objects using different initial 
orbital radii and different giant planet migration rates, 
and concluded that up to 50% of objects could remain in 
the system, albeit with a la rge eccentricity distribution. 
Edga r~fe Artvmowiczl (|2004i ) also examined the dynami- 
cal impact of the migration of a planet of various sizes on 
a planetesimal disk, concluding that eccentricities could 
rise as high as 0.4 but would be damped down quickly in 
the presence of gas. 

Several studies have recently begun to investigate the 
evolution of a more realistic protoplanetary disk, both 



during migration and in a post-migration c o nfigur ation. 
iFogg fc Nelson] (|2005l ) and iFogg fc Nelson! ([20061 ) ana- 
lyzed the survival rate of inner-disk objects during mi- 
gration in a range of size regimes approximating run- 
away growth, including the effects of gas drag and sub- 
sequently a more realistic evolving disk model. They 
concluded that between 50 and 90% of material would 
be retained, with the distribution interior and exterior 
to the giant planet var y ing w ith disk maturity. Alter- 
nately, iRavmond et al.l (|2005f) examined the formation 
and composition of terrestrial planets in the presence of 
an existing inner giant planet at varying stopping dis- 
tances and without a gas disk, concluding that a hot 
Jupiter would have little influence on terrestrial planet 
formation outside its orbit. That study proposed that 
terrestrial planet formation would be inhibited for orbits 
with periods within a factor of 3 - 4 of either an outer or 
inner giant planet. 

In this study we present detailed analysis of numerical 
simulations in which we explore the final stages of ter- 
restrial planet formation after the migration of a giant 
planet. Simulations begin with a two-phase protoplan- 
etary disk and a fully-formed giant planet in the outer 
disk, and follow the system through the giant planet's 
migration and 200 Myr of additional evolution, suffi- 
cient to examine the final characteristics of planets in 
the terrestrial zone. These results bridge the gap be- 
tween the short-term effects of migration on the p lan- 
etesimal disk demonstrated bv lFogg fc Nelson! (|2005| ) and 
the long-term evolution and final configurations of plan- 
etary systems with cl ose-in giant planets explored by 
IRavmond et al.l (|2005l). A subset of t hese si mulations 
were first presented in IRavmond et a l. (2006b) (referred 
to here as Paper I), and in this paper we present two 
additional simulation sets run with different initial con- 
ditions and a detailed description and analysis of all three 
simulation sets. 

In Section [5] we describe the simulation details and ini- 
tial conditions; in Section [3] we describe the results for 
all the simulations and compare the three different mod- 
els; in Section [4] we explore the potential ramifications 
for known extrasolar planetary systems; and in Section [5] 
we provide a summary of our results and suggest future 
work. 

2. MODEL PARAMETERS 

2.1. Initial Conditions 

We performed three different sets of simulations to ex- 
amine the effects of the most important parameters in 
the simulations: the number and location of giant planets 
and the presence of a gaseous disk during and after mi- 
gration. The primary set includes only one Jupiter-mass 
planet that migrates, and includes the viscous damping 
effects of a gaseous disk. The second set has both a mi- 
grating Jupiter-mass planet and a Saturn-mass planet 
stationary at 9.5 AU, with viscous damping included. 
The final set includes the two giant planets but does not 
incorporate gas drag. Our three sets of simulations are 
named JD for 'Jupiter and (gas) Drag', JSD for 'Jupiter, 
Saturn and Drag', and JSN for 'Jupiter, Saturn, No 
drag'. Using different random number seeds, we gen- 
erated five disks of randomly distributed embryos and 
planetesimals for each set. However, one run in each sim- 
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ulation set was unusable due to corruption in the data 
storage or transfers between processors; therefore final 
results are presented here for four simulations in each 
set, resulting in a total of twelve simulation results. The 
first four simulations (JD) were presented in Paper I; we 
incorporate them here as part of the larger set of simu- 
lations. 

The disk models used in this study have three major 
components: a solid disk of small protoplanetary bod- 
ies, a gaseous disk, and one or two fully-formed giant 
planets. We start our simulations at the end of the oli- 
garchic growth phase of the protoplanetary disk, such 
that ^1000-km planetary embryos have formed through- 
out the disk but the mass in km-size d planetesimals 
is compara ble to the mass in embryos (Ko kubo fc Idal 
1998, 2000). Additionally, we incorporate the dissipa- 
tive damping f rom a decaying gas d i sk, which disappear s 
after 10 Myr ([Bricerio et alj 120011: lHaisch et all l2001h . 
Finally, we include either one or two giant planets, as- 
suming them to be fully-formed by this stage due to rapid 
formation processes in the outer disk. This is consistent 
with recent results showing that giant planets may form 
in < 1 Myr via either core-accretion ([Rice fc Armitagel 



I2003t lAlibert et al.ll2005t iHubickvi et all 



tational collapse (Boss 2001]; iMaver et al 



2005) or gra vi- 
12002 . l200l . 



The radial surface density profile of the solid material 
in the disk follo ws from the min i mum-mass sola r nebu la 
model (MMSN: IWeidenschillind (fl97l : lHavashil (fl98lh ): 



(1) 



where Si is the surface density at 1 AU in an MMSN 
disk (~ 6g/cm 2 and 1700 q /cm 2 for t he sol id and gaseous 
components, respectively; lHavashil (|1981h ) and / is a 
scaling factor for the total surface density. Since the 
presence of clos e-in giant planets correlates with stellar 
metallicity (e.g. lLaws et all (|2003[ )), we model a signifi- 
cantly more massive disk than the MMSN, with / as 2.2. 
Similarly, the nebular gas density follows an exponential 
profile of the form 

p(r, z) = p Q (r)exp{~z 2 /z (r) 2 } q/cm 3 (2) 

as described bv lThommes et all ([2003D . where po is the 
midplane density taken from the MMSN: 

p™ in (r) = 1.4 x 10~ 9 (r/l AU)~ 11/4 g/cm 3 (3) 
and zq is disk vertical scale height given by 



Zo(r) = 0.0472(r/l AU) 5/4 AU 



(4) 



The gas density decreases linearly over 10 Myr, sim- 
ulating the removal of the gaseous nebula on ob served 
timescales ([Briceiio et al.ll200U lHaisch et al.ll2001l) . 

Figure [1] illustrates our initial conditions for a single 
run. The initial solid disk is composed of ~ 80 planetary 
embryos and 1200 planetesimals. It extends from 0.25 to 
4.5 AU (with a ~ 2 Hill radii gap on either side of the 
Jupiter-mass planet at 5.2 AU), and then from 6 AU to 9 
AU. The solid disk comprises 17 M®, of which 10 M® is 
equally distributed between planetesimals and embryos 
from 0.25 to 4.5 AU. Embryos between 0.25 and 4.5 AU 
are spaced randomly by A=5-10 mutual Hill radii, and 
have masses between roughly 0.01 and 0.4 M® with a 



mean of 0.11 M®. Embryo masses M em b increase with 
radia l distance r as M em b cx A 3 / 2 r 3 / 4 (j Raymond et al.l 
l200l . In the inner region between 0.25 and 4.5 AU are 
1000 planetesimals of 0.005 M® each, distributed radi- 
ally as r~ 1 ' 2 , i.e. as the annular mass in our surface den- 
sity distribution. The surface density distribution has a 
jump immediately past the snow line, assumed to lie at 
5 AU; therefore, isolation masses are larger in the region 
from 6 to 9 AU. In this outer region we space embryos 
by A=3-6 mutual Hill radii, forming only 4-7 embryos 
between 0.2 and 1.1 M®, with a mean of 0.6 M® and to- 
taling 5.7 M®. In addition, 200 planetesimals comprising 
1.3 M® are placed in the region. The higher embryo mass 
- planetesimal mass ratio in the outer region reflects the 
faster embryo formation time in this higher density re- 
gion. Indeed, massive embryos must have formed quickly 
in this scenario (if giant planets form via core-accretion), 
because the giant planets are already fully-formed. Start- 
ing eccentricities are randomly selected up to 0.02 and 
inclinations are set at 0.1°, but these initial values are 
inconsequential since the distributions are quickly per- 
turbed. 

Each embryo and planetesimal is assigned a composi- 
tion based on its starting location. Water and iron con- 
tents are based on values from our Solar System, where 
comets are thought to contain almost half their mass in 
water ice and asteroids bey ond ~ 2.5 AU co ntain sig- 
nificant quantities of wate r (I Abe et all (120001): se e Sec- 
tion O and Fig. 2 from iRavmond et all (|2004l )). In- 
side 2 AU, embryos and protoplanets are assumed to be 
dry, and beyond 5 AU they contain 50% of their mass 
in water. Between 5 AU and 2.5 AU they contain 5% 
water by mass, and from 2-2.5 AU they contain 0.1% 
water by mass; this distribution corresponds to starting 
mean water mass fraction of ~ 8 x 10~ 3 inside 5 AU 
(see Figure [T]) . Starting iron contents are interpolated 
between the known valu es of the planets (values from 
lLodders fc Feglevl (|1998f) ). including a dummy value of 
0.4 in place o f Mercury because of i t s ano malously large 
iron content ([Raymond et al.ll2004L l2005h ; the starting 
mean iron mass fraction is ~0.32 inside 4.5 AU and ^0.13 
beyond 6 AU. 

Embryos and giant planets interact gravitationally 
with every body in the system. Planetesimals, however, 
are non self-interacting and are given an effective mass 
to simulate a collectio n of much smaller particles (e.g. 
iThommes etall ((2003)). In this way we can realistically 
include the effects of 1) gas drag on planetesimals and 2) 
dynamical friction of planetesimals on embryos while us- 
ing a reasonable particle number. Collisions are treated 
as inelastic mergers conserving wat er and iron (for a dis- 
cussion, see IRavmond etall ([2004D ). 

2.2. Simulation Details 

Numerical simulations were performed using a mod- 
ified version of the publicly-available hybrid sym- 
plectic integr ator package MERCURY by Chambers 
fc M igliorini ([Chambers fc Migliorim] 119971 : IChambersl 
1999). We integrated each simulation for 200 Myr 
with a 2-day timestep, which accurately integrates or- 
bits with apocenter distances larger than roughly 0.05 
AU |L evison fc Duncarj |2000D . Energy is not properly 
conserved for orbits inside 0.05 AU - these bodies are 
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typically given an artificial energy "kick" that results in 
a close encounter and dynamical ejection. This effect is 
important in certain instances and is discussed further in 
Section [331 

To examine the effects of migration, we modified the 
integrator to accommodate an art ificial inward migration 
of a g iant planet, as discussed in iMandell fc Sigurdssonl 
We use a simple drag force as described in 



ot a g i 
(12001 



Chian g~et alj (|2002f ) that produces a linear inward mi- 
gration: 



-Mpvp 



(5) 



where Mp and vp relate to the planet and t m i g is the 
migration time. This method for simulating migration 
produces no artificial changes in eccentricity and incli- 
nation for the migrating planet, and therefore removes 
the potential for non-physical orbital excitation of the 
giant planet. Migration rates as a function of time were 
modeled after recent simulations of Type II giant planet 
migra tion in the literature, which range from 10 5 to 10 6 
years <|D'Angelo et all 120031 : iNelson fc Papak>izoull2004f) 
with an exponential tail to simulate potential braking 
processes (see Sect[T|). 

In addition, a basic fluid drag force was instituted to 
simulate the dissipation of energy due to the surrounding 
gaseous disk. The gas drag takes the basic form of Stokes 
drag: 



a drao — Kv re i~V 



rel (6) 

where v re i is the relative velocity of the object with 
respect to the surrounding gas. The drag parameter is 
defined as 



K = 



^PgasCl 



(7) 



where p gas refers to the local gas density, p m and r m 
refer to the density and mass of the object and the dra g 
coefficient Cd is defined to be 1 (|Adachi et al.l Il976h . 
The gas disk is assumed to revolve in a circular orbit at 
sub-Keplerian velocity due to internal pressure support, 
following the relation 



gas 



= v K (l - 77) 



with 77 defined as 



16 V 



20 0") 



(8) 



(9) 



where a and (3 are the exponents of radial dependen- 
cies of density and temperature respectively and zq is the 
radial ly-dependent vertical disk scale height, defined be- 
low <|Th ommcs et all 120031 ). The solid bodies therefore 
experience both a orbital damping effect and an inward 
migration, decreasing with mass and increasing with gas 
density. 

In these simulations we do not account for torques on 
sub-Saturnian sized bodies due to density waves in the 
gas disk, which may lead to a reduced velocity disper- 
sion ( known as gravitational or tidal gas drag; iWardl 
(1993) ) and/or or bital decay (known as Type I migra- 
tion: IWardl (|1997f) ) of embryo and planet-sized bodies. 



We neglect these effects primarily because the role of 
disk torques on fully-embedded objects is still uncertain. 
Type I migra tion rates from simulations range from 10 3 
to 10 6 years (jTanaka et al.ll2002t Ir/Angelo et all 120031) . 
which would alternately result in either the loss of all 
solid bodies before the gas disk has disappeared, or very 
little change in orbital position for large bodies. It 
is even unclear whether migration or damping due to 
disk t orques functions at all in realistic turbulent ga s 
disks (INelson fc Papaloizoul 12004 lLaughlin et al.ll2004f) . 
iMcNeil et al.1 (|2005f ) investigated the oligarchic stage of 
terrestrial planet formation for different Type I migra- 
tion scenarios, demonstrating that faster migration rates 
serve to damp down excitation and widen embryo spac- 
ing for material within ~ 2 AU. This may prove to be 
important in determining the final mass and composition 
of terrestrial planets, but since the rate of migration of 
protoplanets compared with the migrating giant planet 
has a significant impact on the scattering energy and 
angular momentum transfer between the giant planets 
and the surrounding solid bodies, arbitrarily choosing a 
migration rate would strongly influence the results. We 
therefore chose not to include it in this work. Future 
simulations will incorporate more detailed models of em- 
bryo migration taken from improved hydrodynamic sim- 
ulations. 

3. RESULTS 

The evolution of the protoplanetary disk can be sep- 
arated into distinct dynamical stages. Migration (Stage 
1) encompasses the migration period, during which the 
giant planet moves through the interior regions of the 
disk. Gas Dissipation (Stage 2) begins after giant-planet 
migration has ceased, and continues until the remaining 
gas disk has fully dissipated. Clearing (Stage 3) includes 
the final stages of accretion and clearing in the absence 
of any gaseous material. In the following subsections 
we describe each stage, and illustrate the differences be- 
tween the results for the three different models used in 
our simulations. Snapshots during the different stages of 
the evolution of each type of simulation are illustrated 
in Figures 1 (JD), 1 (JSD), and[4] (JSN). Each plot fol- 
lows the evolution of all bodies in the simulation, and 
tracks their masses, orbital parameters and compositions 
through time. In addition, the mean and range of orbital 
and compositional properties for each model at the end 
of each stage are listed in Tables [TJ [31 and [31 

3.1. Stage 1: Migration 

During Stage 1 the Jupiter-mass planet migrates in- 
wards from 5.2 to 0.25 AU over 10 5 years, and the inward 
progression of mean motion resonances results in the 
shepherding and orbital excitation of the disk material 
interior and exterior to the giant planet. As the Jupiter- 
mass planet continues through the inner disk, material 
is either scattered outward or captured into resonant or- 
bits and forced to migrate inward with th e Jupiter-mass 
planet . This resonance capture, noted bv iTanaka fc Idal 
(1999), occurs when a body close to a strong mean mo- 
tion resonance (e.g., the 2:1 MMR) has its orbit excited 
and its eccentricity increased. Thus, the body's perihe- 
lion distance has decreased and it has lost angular mo- 
mentum. An eccentric orbit increases the relative veloc- 
ity between the protoplanet and the pressure-supported 
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gas disk, thereby increasing the effects of gas drag. In ad- 
dition, an eccentric orbit encounters more nearby plan- 
etesimals and also enhances dynamical friction. These 
dissipative forces tend to reduce the body's energy and 
decrease its semi-major axis and eccentricity. Thus, the 
protoplanet is moved inward of the resonance with the 
giant planet; as the giant planet continues to migrate, the 
resonance can continue pushing it inward due to these an- 
gular momentum loss mechanisms. Smaller bodies that 
feel gas drag more strongly can be pushed inward by 
higher-order resonances. This is clearly seen in Fig. [2l 
where protoplanets are pushed inward by the 2:1 and 
3:2 resonances, and planetesimals by the 8:1 resonance. 
The decreasing spacing between resonances enhances the 
accretion and scattering rate as the giant planet moves 
inward, with planetesimal encounters peaking at approx- 
imately 7 x 10 4 years when the Jupiter-mass planet's 
semi-major axis reaches ~ 0.4 AU and the 2:1 resonance 
reaches the inner edge of the disk (see Figure [5]). At the 
end of the giant planet's inward migration, the remaining 
disk material is divided between bodies captured in very 
close low eccentricity orbits in interior resonances with 
the Jupiter-mass planet and bodies in high eccentricity 
orbits beyond 0.5 AU. 

In the simulations including the presence of a gas disk 
(Models JD and JSD) orbital excitation of planetesimals 
is almost completely damped in the inner disk - by the 
end of Stage 1 the mean eccentricity for planetesimals 
within 9.5 AU is less than 0.1, compared with approxi- 
mately 0.5 for embryos; mean inclination follows similar 
trends but with slightly higher values for bodies in the 
JSD simulations due to multiple scatterings (see Table 
[1]). The excitation of embryos is damped by a combina- 
tion of relatively weak gas drag and dynamical friction 
from interactions with planetesimals. These effects are 
sufficiently strong to halt almost all ejections, and the 
orbital radii of scattered bodies are limited to within 50 
AU. Additionally, these forces result in strong resonant 
trapping by the migrating planet, frequent collisions and 
a rapid increase in embryo masses - Figure [7] reveals that 
the large majority of both embryo-embryo collisions and 
embryo-planetesimal collisions occur during the migra- 
tion period, and Figure [5] demonstrates the loss of plan- 
etesimal mass to embryo mass. By the time migration 
has ended, the combination of rapid accretion and scat- 
tering has cleared almost all planetesimals from within 
5 AU, and the average embryo mass for this region has 
risen from 0.11 M ffl in ~65 objects to 0.45 Mq in ~17 
objects. Eccentricities for remaining embryos beyond 0.5 
AU are evenly distributed from 0.2 to 0.8, increasing with 
semi-major axis (see Figures [2] and |3]) ; this relationship 
is charac teristic of a scattered populat ion of embryos as 
noted bv lMandell fc Sigurdsson! (|2003l ). Inclinations re- 
main low due to damping, with average inclinations re- 
maining below 5° (see Tabled]). Within 5 AU the mean 
composition of embryos does not change drastically since 
the bodies surviving in the inner disk are not augmented 
by any additional material over these short timescales, 
but the material is locally mixed. The exception is the 
JSD model, where water-rich material is scattered inward 
by the outer Saturn-mass planet and the mean water 
mass fraction increases by a factor of 2. 

Beyond 5 AU the dynamical effects of the presence of 
a second giant planet become important. In simulations 



without a second planet, the only external excitation for 
objects orbiting beyond the initial orbit of the migrat- 
ing Jupiter-mass planet comes from external resonance 
crossings and mutual gravitational interactions. Material 
native to this region remains relatively undisturbed, with 
very little orbital excitation. Bodies scattered into the re- 
gion by the migrating giant have large eccentricities that 
are subsequently damped on Myr timescales. Maximum 
orbital eccentricity for these outer embryos increases up 
to ~ 0.5, and by the end of migration dynamical fric- 
tion and damping due to gas drag reduce the maximum 
eccentricity to ~ 0.2 (see Figure [2|). In the JSD simu- 
lations, small bodies in the outer system are excited by 
both mean motion resonances and gravitational scatter- 
ing. This results in a significant scattered planetesimal 
population beyond 9.5 AU with a mean eccentricity of 0.3 
and inclinations up to 30°. A small number of embryos 
are also scattered by both the migrating Jupiter-mass 
planet and the outer Saturn-mass planet, resulting in a 
semi-major axes as high as 50 AU. 

In addition to the presence of material remaining out- 
side the orbit of the close-in giant planet, all simula- 
tions including gas drag form one or more terrestrial- 
mass bodies residing within the orbit of the Jupiter-mass 
planet (see Figure [T0|) . These small bodies are the result 
of resonant shepherding by the giant planet as it migrates 
inward, and accretion at the resonance can build up bod- 
ies of several M®. These "Hot Earths" are analogous 
to the 7.5 planet found by radial velocity searches 
around the M star GJ 876 (Rivera et al. 2005) and sug- 
gest that this phenomenon is common. Other computa- 
tional studies have demonstrated reso nant shepherding 
as a mechanism to explain hot Earth s (|Zhou et afll2005l : 
iFogg fc Nelsonl 120051 . I2006T) . though iZhou et all (l2005h 
explain the planet around GJ 876 by invoking moving 
secular resonances during disk dissipation rather than 
mean-motion resonances during migration. However, in 
almost all of the simulations the inner 'hot Earth' be- 
comes unstable as the gas disk disappears and either im- 
pacts the Jupiter-mass planet or is ejected. Similarly, 
in simulations without any gas drag almost no accretion 
occurs interior to the migrating giant pl anet. This sen- 
sitivit y to damping forces was noted bv IFogg fc Nelsonl 
(2006), but may be mitigated by the presence of colli- 
sional debris not included in these simulations. More 
detailed work is necessary to accurately model the long 
term stability of these bodies. 

Simulations performed without the presence of a gas 
disk (Model JSN) demonstrated very different encoun- 
ters between bodies than Models JD and JSD through- 
out the disk. Without the damping effects of gas drag 
on planetesimals, accretion ceases after ~ 40, 000 years 
as the Jupiter-mass planet enters the inner system (see 
Figure [7|) . Excitation from resonance crossings and di- 
rect scattering cause ~ 40% of the mass in small bodies 
to be lost from the system by the end of migration (see 
Figure [5] and Figure [6]) , with the remainder of material 
distributed evenly in eccentricity and inclination phase 
space. Approximately 65% of the remaining mass still 
resides within 10 AU (compared to 75% for simulations 
with gas drag), but the ratio of the total mass in embryos 
to the total mass in planetesimals in this region has only 
risen to only ~ 2, compared to ~ 10 for simulations in- 
cluding the effects of a gaseous disk. The mean embryo 
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mass is 0.18 Mq compared with 0.45 M® for the sim- 
ulations with gas drag, and the mean eccentricity and 
inclination for both embryos and planetesimals are both 
significantly higher than the JD and JSD simulations. 

3.2. Stage 2: Gas Dissipation 

Once the Jupiter-mass planet has ceased migrating at 
0.25 AU, the remaining bodies in the disk continue to 
scatter and accrete each other during Stages 2 and 3, 
similar to the classical oligarchic and chaotic stages of 
terrestrial planet formation. At the start of Stage 2 the 
protoplanetary embryos in the disk are dynamically hot, 
with mean eccentricities increasing with orbital radius. 
However, during Stage 2 the presence of gas drag and 
dynamical friction serve to damp down eccentricities and 
inclinations and replenish the inner disk due to migration 
of planetesimals. During this stage the location and rate 
of accretion depends critically on the balance between 
the size of planetesimals and embryos and the density of 
the gas disk: faster migration leads to higher accretion 
rates at smaller orbital radii, while damping of dynam- 
ical excitation leads to reduced accretion and resonance 
capture. 

After a hiatus of several million years due to clearing by 
the giant planet migration and damping by the remain- 
ing gas disk, the accretion of planetesimals increases over 
the final duration of Stage 2 as a result of the refilling 
of feeding zones due to planetesimal orbital decay and 
increased orbital crossing due to decreased gas drag (see 
Figure [7] for the accretion history of the combined data 
sets). The accretion wavefront proceeds inward as the 
gas becomes less dense and orbital excitation increases. 
This late-stage accretion is responsible for the high water 
mass fractions seen in the final configurations - embryos 
acquire most of their mass during Stage 1, but most of 
their water during Stage 2 (see Figure [9]). In simulations 
with the JD model accretion continues until after the 
gas has disappeared, but the presence of a second mas- 
sive planet in the JSD model leads to clearing between 5 
- 15 AU and a decrease in accretion after 5 x 10 6 years. 
The outer planet also serves to clear out high-eccentricity 
embryos, leading to a much lower mean eccentricity and 
fewer embryos in the outer disk as compared with the JD 
model (see Table [2]). After 10 Myr the mass in planetes- 
imals in both the JD and JSD models is only ~ 1M®, 
but the total mass in the JSD simulations is less than 
30% of the mass remaining in the JD simulations due to 
embryo ejection in the outer system (see Figure [SJ Table 
[2]). Eccentricities and inclinations in the inner system 
for both models decrease dramatically over the lifetime 
of the gas disk, with embryo eccentricities and inclina- 
tions decreasing by 20-50% and planetesimals attaining 
a mean eccentricity below 0.1 and a mean inclination less 
than 2°. 

3.3. Stage 3: Clearing 

In the simulations with gas drag, Stage 3 begins as the 
final gas disk dissipates and dynamical interactions be- 
gin to excite the remaining solid bodies. Accretion rates 
decrease steadily and cease by approximately 30 Myr, 
and scattering and ejection clear out most of the mate- 
rial in the inner disk over the remaining 170 Myr (see 
Figure [6]). At the end of the simulations only a few large 



protoplanets are left in the inner disk, but longer dynam- 
ical timescales in the outer disk result in longer clearing 
times and increased final eccentricities and inclinations 
for the JD simulations - the mean eccentricity rises to 
~0.3 and inclinations increase up to ~50° (see Table 
In contrast, the Saturn-mass planet begins to excite and 
scatter bodies near the end of Stage 2, and continues to 
eliminate highly excited bodies so that the mean eccen- 
tricities and inclinations for embryos at the end of 200 
Myr are below 0.2 and 20° respectively. The few remain- 
ing planetesimals are scattered randomly in eccentricity 
and inclination space. 

In the JSN simulations, the system moves directly from 
Stage 1 to Stage 3. Ejections of both embryos and plan- 
etesimals increase over the first 10 Myr, then begin to tail 
off (see Figure H]). Eccentricities maintain a flat distribu- 
tion up to 0.7 due to a random distribution of scattering 
energies, while inclinations are spread out to 60°. More 
than 70% of the total mass is lost by 5 Myr due to ejection 
and impacts with the central star, leaving only a small 
number of embryos in chaotic orbits and a scattered pop- 
ulation of planetesimals (see Figure [4]). By 10 Myr the 
total mass in embryos is only 3 times the mass in plan- 
etesimals, compared with a total embryo/planetesimal 
mass ratio greater than 10 for the simulations with gas 
drag included. Highly excited bodies are cleared out of 
the inner system over the remainder of the simulation 
period and the mean eccentricity and inclination of the 
remaining embryos decrease to 0.25 and 17° respectively, 
but the planetesimals remain in highly inclined and ec- 
centric orbits. 

3.4. Final Configurations 

Figure QT] shows the final configurations of all twelve 
of our simulations, with the Solar System included for 
comparison. The details of the final planetary properties 
for bodies within 5 AU are listed in Tables |H E ] and E H 
As seen in prev i ous sim ulations (e.g.. IChambersI (|2001[ ): 
I Raymond et al.l (|2004f )). a diversity of outcomes exists 
within each group of similar simulations. 

Clear differences can be seen between the different sets 
of simulations, in particular between simulations with 
(JD and JSD) and without (JSN) gas drag. Both the 
mean planet mass and the mean water mass fraction in 
simulations with drag is much higher than without drag. 
These trends have two root causes: 1) the survival of 
bodies scattered during migration depends strongly on 
drag as a dissipative force, leaving more material in the 
system when gas is present and promoting radial mix- 
ing, and 2) icy planetesimals from beyond 5 AU can spi- 
ral in to the terrestrial zone over the lifetime of the gas 
disk. Simulations with drag also form large hot Earths 
due to resonance capture by orbital damping (discussed 
in Section 13. 1|) ; Figure [10] demonstrates that all the hot 
Earths lie just inside resonances with the hot Jupiter. 
These hot Earths are rare in simulations without drag - 
only simulation JSN-3 formed a dry, 0.22 M® hot Earth 
at 0.05 AU. We attribute the existence of this planet to 
damping from dynamical friction during migration rather 
than gas drag, but scattering and ejection is much more 
likely under these circumstances and the probability of a 
massive body surviving is low. Also, the final radius of 
the close-in giant planets' orbits are scattered in the gas 
drag simulations, due to the artificial migration caused 
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by non-physical giant planet - planetesimal effects (dis- 
cussed in Section 13. 5j) ; in simulations without drag the 
giant planet remains at 0.25 AU. 

There also exist differences between simulations with 
and without an outer giant planet. The JD simula- 
tions contain extensive scattered disks beyond a few AU 
with high inclinations and eccentricities and very long 
accretion timescales, reminiscent of the Sola r System's 
scatt e red disk of comets bey ond Neptune (|Luu et al.1 
fl997t |D uncan fc Levisonl |1997( ). These disks are likely 
to produce vast amounts of cold dust and may be de- 
tectable around other stars; late-stage accretion and 
fragmentation due to this scenario could be responsible 
for the debris disks seen around intermediate-age main- 
sequence stars such as (3 Pictoris (Golimo wski et aTll2006b 
iBrvden et al. 2006) . Outer disk material is cleared out in 
simulations with an outer giant planet, greatly reducing 
the amount of surviving outer disk material, suggesting 
that dust mass production could correlate positively with 
giant planet migration and negatively with giant planet 
survival in the outer system. 

Potentially habitable planets survive in almost all of 
the simulations. In the eight simulations with drag, seven 
planets larger than 0.2 Mq (including two in simulation 
JSD-4) formed in the Habitable Zone, the orbital region 
where the stellar flux is sufficient to maintain liquid wa- 
ter on the surface of a planet (ass umed to lie between .8 
and 1.5 AU for these simulations; iKasting et alj |l993)). 
Five planets with masses from 0.13 to 0.4 Mq formed in 
the Habitable Zone in the four simulations without drag. 
In two separate cases, however, two planets remain in 
the Habitable Zone on crossing orbits (sims. JSN-2 and 
JSN-5). In each case both of the planets have inclina- 
tions of ~ 20° so the timescale for the two planets to col- 
lide is long. Long-term dynamical studies are necessary 
to confirm the stability of these planets on billion-year 
timescales. Planets as large as 1 M ffi form within 5 AU 
in simulations without gas drag, but it seems clear that 
additional material (either through a more massive disk 
or an inward flux) is necessary to form similarly sized 
planets in the Habitable Zone. 

3.5. Drag- Induced Inward Migration 

In some cases, the inner gas giant and/or surviving em- 
bryos underwent additional migration during the course 
of Stage 2. The cause of this additional migration origi- 
nates with the resonance capture of planetesimals spiral- 
ing inward due to gas drag. In most cases, the planetes- 
imals become trapped in outer resonances, then cascade 
down to just outside the 2:1 mean motion resonance with 
the giant pla net. This o uter r esonance capture was in- 
vestigated by lThommesI (|2005| ) with regard to larger ob- 
jects, focusing on differential migration due to differences 
in Type II (giant planet) and Type I (terrestrial planet) 
migration rates. In this case, the inward migration comes 
from aerodynamic drag, and the capture process is aided 
by damping of eccentricities by the gas disk. The ef- 
fect is strongest in the inner disk where gas densities are 
high, and once the gas disk dissipates the planetesimals 
are excited and either accreted onto the larger object or 
ejected from the system. 

The combination of angular momentum loss from the 
larger object to the planetesimals (through orbital exci- 
tation from the strong resonance interaction) and trans- 



fer of this angular momentum from the planetesimals to 
the gas disk (in the form of orbital damping) leads to a 
net loss of angular momentum from the solid bodies and 
subsequent orbital decay. The amount of inward migra- 
tion varied from simulation to simulation, depending on 
the orbital radius and number of planetesimals trapped 
in resonance. The total amount of inward movement of 
giant planets ranged from zero to almost 0.2 AU - in 
most cases it was about 0.1 AU. For terrestrial embryos 
at 1 AU the inward migration was similar, but no ad- 
ditional inward motion was apparent for objects beyond 
1.5 AU (see Figure [8] for examples). The effect on the 
final orbits of the terrestrial bodies was usually overshad- 
owed by orbital rearrangement due to late impacts and 
scattering events. The additional inward motion of the 
inner Jupiter-mass planet also caused the orbits of close- 
in terrestrial planets in some simulations to become too 
small to be properly reso lved with our 2 day timestep 
(jLevison &; Duncan! (|2000f )). These orbits became unsta- 
ble and usually collided with the close-in giant planet. 

Though this additional source of inward migration rep- 
resents an interesting avenue for future work, the effect 
cannot be properly analyzed in this study due to our un- 
usual treatment of the physical qualities of the planetes- 
imals. As noted in Section 12. 2[ the planetesimals have 
a dynamical mass of a body with a size of ^1000 km, 
but interact with the gas disk with a more realistic cross 
section of a collection of 10 km bodies. Thus, the rate 
of angular momentum transfer for a single planetesimal 
caught in resonance is equivalent to that of 10 6 smaller 
bodies. Additionally, for computational efficiency plan- 
etesimals do not self-interact. Thus, when trapped in 
the same resonant orbit these small bodies are unable 
to coagulate into a larger body that would lose angular 
momentum to the gas disk less efficiently. The angular 
momentum loss and subsequent orbital decay over the 
duration of the gas disk is therefore much greater than 
would be present in a realistic system. 

3.6. Water Delivery 

Water-rich material is accumulated by the growing 
planets throughout their accretion. As shown in Paper 
1, two primary mechanisms contribute: 1) radial mix- 
ing induced by large eccentricities as a result of scat- 
tering interactions with the migrating giant planet; and 
2) accretion of in-spiraling icy material because of gas 
drag. The first mechanism delivers water mainly in the 
form of hydrated asteroidal material and takes place 
through each stage, although the vast majority of ra- 
dial mixing occurs through giant planet scattering dur- 
ing Stage 1, with additional embryo-planetesimal scat- 
tering occurring during Stage 3 after the gas disk has 
dissipate d. This mode of tran sport is similar to the 
model of iMorbidelli et ail (|2000j ), except that it is sub- 
stantially enhanced by the large eccentricities induced by 
close encounters between the migrating giant planet and 
the smaller bodies. The second mechanism occurs mostly 
during Stage 2 as water-rich material filters through the 
inner system over long timescales. In our own system this 
mechanism most likely contributed little if any water to 
the Earth because of Jupiter's presence: in-spiraling icy 
bodies would have first encountered Jupiter and been 
scattered outward or accreted. In the simulations pre- 
sented here, the very water-rich planetesimals originat- 
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ing past 5 AU move inward quickly, but damping of ex- 
citation diminishes collisions with embryos and allows 
the small bodies to reach sub-AU radii. Accretion of 
these bodies tends to occur shortly after gas dissipation 
when damping ceases and their eccentricities and incli- 
nations are increased by perturbation from nearby ter- 
restrial bodies, resulting in collisional orbits. 

As discussed in Paper 1 and below, the water contents 
of Habitable Zone planets in the simulations with gas 
drag are very high, with typical values of ~ 10% water 
by mass. However, it is difficult to translate these into 
likely va lues, since we do not take water depletion fr om 
impacts (jGenda fe Abdl200alCanup fc Pierazzol2006D or 
hydrodynamic escape ( Matsui fc Abdl 19861 ) into account. 
Previous simulations with a single, non-migrating outer 
giant planet form planets in the Habitable Zone with 
~5x 10 -3 water by mass (Raymond et al 2004, 2006a). 
If we assume that previous simulations would have re- 
sulted in water contents similar to the Earth's (~ 10~ 3 
by mass) if depletion were accounted for, then the present 
simulations should still have about 20 times as much 
water as the Earth. However, it is possible that previ- 
ous simulations instead formed planets which would have 
more water than the Earth because the giant planets 
in those simulations were on circular orbits rather than 
mildly eccentri c orbits similar to the orb i ts of our own gi- 
ant p lanets (|Chambers fc Cassenl l2002t iRavmond et al.l 
2004). Thus, our current planets may in fact contain 
significantly more water and are likely to be covered in 
global oceans several km deep and be veritable "ocean 
planets" (e.g. lLeger et al.l (|2004f )). 

The water contents of Habitable Zone planets in simu- 
lations with and without gas drag differ by about two or- 
ders of magnitude. The reasons for this contrast are tied 
to the presence of nebular gas, via radial mixing of sur- 
viving bodies and in-spiraling of icy planetesimals. Thus, 
the water contents of our final planetary systems may in 
fact be dependent on our assumptions about the lifetime 
and depletion of the gas disk and the effective mass of 
planetesimals. For instance, if gas disk dissipation is the 
halting mecha nism of the giant pla. net migration (as as- 
sumed by, e.g.. lTrilling et all (|1998l )1. then perhaps only 
a Jupiter-mass or so of gas remains at the end of migra- 
tion. We expect that re-circularization of scattered mate- 
rial and inward migration of small bodies would continue 
to take place to some extent in such a disk, such that a 
significant fraction of material would survive to repopu- 
late the terrestrial zone. Thus, the amount of material 
accreted from the outer system by planets in the inte- 
rior of the disk would be less. The terrestrial planets 
in such systems would probably represent an intermedi- 
ate case between our simulations with and without drag, 
with moderate masses and water contents. As shown by 
the final configurations, we can expect that such systems 
would form habitable planets similar to those formed 
here. 

4. EARTH-LIKE PLANETS IN KNOWN 
PLANETARY SYSTEMS 

It has been shown that an Earth-like planet would be 
dynamically stable for long timcscalcs in the Habitable 
Zones of many of the known sy st ems of giant plane ts 
(e.g., iMenou fc Tabachnikl (|2003f) ; Uones et alj (|2005l) ). 
However, these studies only test the stability of existing 



particles over long timescales; it is assumed that initially- 
stable objects would form in-situ. This is not the case in 
moderately unstable environments; stability models pre- 
dict the existence of an Earth-mass planet at ~ 3 AU 
in the our system's Asteroid Belt, where such a planet 
would indeed be dynamically stable. Clearly, when as- 
sessing the potential for a viable location for habitable 
planets it is critical to model the actual formation pro- 
cess. 

Our results can be applied to constrain which of the 
known systems of extra-solar giant planets could have 
formed a terrestrial planet in the Habitable Zone. By 
combining limits on formation of terrestrial planets in 
the presence of various giant planet scenarios, we can 
make a more rigorous assessment of the possibility for 
terrestrial planets in known planetary systems and begin 
to develop criteria for a habitable system. In this section 
we combine limits on the formation of terrestrial planets 
in the presence of various giant planet scenarios from this 
study and a previous study, and assess of the probability 
of the formation of terrestrial planets in known planetary 
systems. 

Examining the typical spacing between terrestrial and 
inner giant planets in our simulations, we find that a 
giant planet must be interior to approximately 0.5 AU 
for a planet at t he low er limit of habitability (~ 0.3 M®; 
iWilliams et al.l (|1997f) ) to form inside the outer bound- 
ary of the Habitable Zone (we have assumed the outer 
edge of the Habitable Zone to be 1.5 AU, although 
this value is uncertain and depe nds on the conden- 
sation properties of CO2 clouds dKasting et al.l [19931 : 
iForget fc Pierrehumbert] 119971 ; iMischna et al.l 12000)). 
The spacing of planets that form exterior to the close- 
in giant varies a great deal from simulation to simula- 
tion. In dynamical terms, what is relevant is the ratio of 
orbital periods of the innermost planet with significant 
mass (>0.3 M$) and the close-in giant planet. This 
value ranges from 3.3 to 43 in our simulations, with a 
mean [median] of 12 [9]. Although there clearly exists a 
range of outcomes, we can define a rough limit on the or- 
bit of an inner giant planet that allows a terrestrial planet 
to form in the Habitable Zone (the most optimistic case, 
i.e. the closest spacing, puts the giant planet at about 
0.7 AU). 

We suggest that terrestrial planets of significant mass 
can form in the Habitable Zone of a Sun-like star if no 
giant planets exist between 0.5 and 2.5 AU. We derive 
these values using the reasoning above based o n our cur- 
rent r esults, in combination with the results of lRavmondl 
(2006), who used hundreds of low-resolution accretion 
simulation to constrain which outer giant planets can 
form Habitable Zone planets (see Fig. 2 from that pa- 
per). These values are clearly eccentricity dependent - 
the characteristic spacing between giant and terrestrial 
plane ts increases quic kly with giant planet eccentricity 
(e.g.. IRavmond! (|2006| )). The highest eccentricity of the 
close-in giant planets in our simulation is about 0.1, and 
this value is similar for most simulations. For our inner 
limit, we therefore assume that a system with a giant 
planet inside 0.5 AU and with an eccentricity less than 
0.1 can form a planet in the Habitable Zone. 

We apply these rough limits to the known extraso- 
lar planets. Our sample consists of 207 planets in 178 
planetary systems (with 21 multiple planet systems) dis- 
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covered by Aug 1 2006. We include planets discovered 
via the radial velocity, micro-le nsing, transit and d irect 
imaging techniques (data from IButler et alj (|2006l ) via 
exoplanets.org, and ISchneiderl (|2006l ) via exoplanet.eu). 
We exclude planets more massive than 15 Jupiter masses 
(Mj; roughly the brown dwarf limit), unless they are 
part of multiple systems. Because the extrasolar planet 
host stars have a range in masses, we must calibrate our 
limits. To do so, we use a mass-lumino s ity re lation that 
is a fit to data of iHillenbrand fe Whitel (|2004h : 

y = 4.101a: 3 + 8.162a; 2 + 7.108a; + 0.065, (10) 

where y = logio (L*/L Q ) and x = logio (M*/M Q ) (John 
Scalo 2006, personal communication). We assume a Hab- 
itable Zone of 0.8-1.5 AU around a solar-mass star, and 
assume that its inner and outer limits scale with the stel- 

1/2 

lar flux, i.e., as LJ . We then assume that the dynamical 
scaling between the Habitable Zone and our giant planet 
limits, measured in terms of orbital period ratios, is in- 
dependent of stellar mass. Table [7] lists our inner and 
outer giant planet limits for a range of stellar masses. 

Figure [T^] shows the known extrasolar planets that 
meet our criteria and are good candidates for having a 
terrestrial planet in the Habitable Zone. Potentially hab- 
itable systems are also listed in Table sorted by host 
star mass. Of the 178 systems in our sample, 65 (37%) 
may have formed a terrestrial planet of 0.3 M ffi or more 
in the Habitable Zone and so are considered to be poten- 
tially habitable systems. Seventeen systems satisfy our 
outer giant planet limit; these are Solar System- like in 
that they have a giant planet exterior to the Habitable 
Zone. In contrast, fifty systems are decidedly different 
than the Solar system, with the potential for an Earth- 
like planet coexisting with a close-in giant planet. One 
system, 55 Cancri, falls into both categories, with three 
close-in giants and a distant outer giant. Nine of the mul- 
tiple planet systems contained a planet that satisfied our 
limits but another planet that did not - 55 Cnc was the 
only multiple system for which all planets met our limits. 
If we consider the ensemble of planets as a whole, then 75 
out of 207 planets (36 %) met our criteria. In the case of 
planets at very large orbital dist ances, such as t he on es 
detected by direct imaging (e.g- IChauvin et alJ (|2005T) ). 
the giant planet is so far away that it will be completely 
dynamically detached from the terrestrial region. 

Our estimate that about one third of the known gi- 
ant planet systems can harbor potenti ally habitable plan - 
ets diverges from previous estimates. I Jones etldl (2005) 
used the stability of Earth-mass test particles to sug- 
gest that about half of the known systems could contain 
a planet in the Habitable Zone. As noted above, the 
stability of a given region does not always imply that a 
planet exists in the region - the Asteroid Belt is a conve- 
nient example of a stable region that contains no Earth- 
mass planet. We susp ect that the d ifferen ce between our 
estimate and that of I Jones et al.1 (|2005f l are a number 
of systems in which stable regions exist but planets are 
unlikely to form. This is not unexpected; for example, 
I Barnes fe Raymond! (|2004h found stable regions for test 
particles in four known e xtrasolar planetary systems, but 
I Raymond et al.l (|2006af ) showed that sizable terrestrial 
planets could only accrete in two of the systems. 

Conversely, in previous work studying the formation of 



planets in systems with one giant p lanet exterior to the 
terrestrial planets (Raymond 2006) we concluded that 
only 5% of the known giant planet systems could form 
planets of at least 0.3 M ffi in the Habitable Zone. Our 
current estimate has increased to roughly one third, be- 
cause of two fa ctors: 1) our sample is larger than in 
iRavmondl (|2006f ). and 2) we have determined that many 
inner giant planet systems could form Earth-like planets 
in the Habitable Zone. Given that observational biases 
make it much easier to find close-in giant planets, it is not 
surprising that our estimate for the number of habitable 
systems has increased. 

In this analysis, we have only considered a few key pa- 
rameters - stellar mass, planet semi-major axis and ec- 
centricity. Table [8] also lists the other parameters which 
we consider most important in terms of habitable planet 
formation: planetary mass and stellar metallicity. Per- 
turbations from lower (higher)-mass planets are corre- 
spondingly weaker (stronger), so for less (more) massive 
giant planets our inner limit would move out (in), and 
our outer limit would move in (out). With the increas- 
ing number o f lower-mass e xtraso l ar planets b e ing d is- 
cover ed (e.g., IButler et all (|200l : lLovis et all |2006)), 
we suspect that our limits are relatively conservative for 
the population of extrasolar planets as a whole. The me- 
dian mass of our sample is 0.76 Mj up , so our limits are 
appropriate on a statistical level. The stellar metallicity 
is thought to be related to the gas- to dust ratio of the 
disk and therefore the total mass in terrestrial building 
blocks. In addition, if the disk mass scales with stellar 
mass (as is generally thought), there may be a deficit 
of rocky material and perhaps fewer Earth-lik e planets 
around low- mass stars ((Raymond et al.|[2006df ). 

The next generations of space missions plan to discover 
and characterize Earth-like planets around other stars. 
These missions include NASA's Kepler, SIM and Ter- 
restrial Planet Finder, and CNES's COROT and ESA's 
Darwin. Our results suggest that terrestrial planets can 
coexist with both close-in giant planets and giant planets 
in outer orbits, expanding the range of planetary systems 
that should be searched with these upcoming missions. 
However, transit searches are very sensitive to orbital in- 
clination. In many cases our simulated Habitable Zone 
planets have mutual inclinations of 5-10° with respect to 
the orbit of the close-in giant planet, making detection 
of transits for both planets impossible. Final inclinations 
are smaller for systems with an exterior giant planet, 
making systems with both an exterior and interior giant 
planet more amenable to detection, but overall the prob- 
ability of seeing the transit of a close-in giant planet and a 
HZ planet in the same system is small. However, transit 
timing measurements may be able to infer the presence 
of HZ planets in these systems (|Holman fc Murravll2005l : 
lAeol et al.lf2005h . 

5. CONCLUSION 

5.1. Summary 

Our simulations demonstrate that terrestrial accretion 
can occur during and after giant planet migration on 
several fronts: 

1. Interior to the migrating giant planet, planetesi- 
mals and embryos are shepherded inward by the 
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combined effects of moving mean motion reso- 
nance s and dissipation via gas drag (iTanaka fc Idal 
(11999ft: also show n by IFogg fc NelsonI (|2005ft and 
IZhou et all ((2005)). Rapid accretion occurs over 
the duration of migration (Stage 1), and 40% — 50% 
of material inside 5 AU typically ends up in the 
form of 1 - 3 "hot Earths" when a gaseous disk 
is present. These planets have masses up to 5 M ffi 
and accrete on the migration timescale of 10 5 years. 
Hot Earths tend to lie inside strong resonances 
such as the 2:1 and 3:2. Planetesimals, which feel 
stronger drag, can be shepherded by higher-order 
resonances. 

2. Exterior to the giant planet, scattered embryos and 
planetesimals have their orbits re-circularized by 
damping from gas drag (felt more strongly by plan- 
etesimals) and dynamical friction (felt by embryos 
due to the planetesimals) during Stage 2. The em- 
bryos in this region begin with compositions simi- 
lar to the initial mean composition, but an influx 
of material into the terrestrial zone (between the 
close-in giant and 2.5 AU) from orbital decay of 
planetesimals due to gas drag can significantly en- 
hance the masses and fraction of water-rich ma- 
terial in these bodies. From these protoplanets a 
system of terrestrial planets forms on a ~ 10 s year 
timescale after the gas disk has dissipated (Stage 
3). 

3. In the outer disk, damping via both gas drag and 
dynamical friction is much weaker. Thus, a scat- 
tered disk of planetesimals and embryos remain 
on high eccentricity and inclination orbits, sim- 
ilar to the scattered disk of comets in the sola r 
system (|Luu et alj|1997t iDuncan fc Levisonlll997ft . 
Timescales for accretion are very long. The pres- 
ence of a second giant planet serves to remove dy- 
namically hot bodies, diminishing the orbital exci- 
tation of bodies and reducing the clearing time in 
the outer system. 

Planets formed in the Habitable Zone in a signifi- 
cant fraction of our simulations. These planets have 
masses and orbits similar to those seen in pr evious sim- 
ulations including only ou t er gian t planets (|Chambersl 
I200H iRavmond et al.ll2004 l2006d) . However, Habitable 
Zone planets in systems with close-in giant planets tend 
to accrete a much larger amount of wat er than those in 
system s with only outer giant planets (|Ravmond et al.l 
l2006bft . The reason for these high water contents is 
twofold: 1) strong radial mixing is induced by the giant 
planet's migration, and 2) in-spiraling icy planetesimals 
are easily accreted by planets in the terrestrial zone. Al- 
though we have not taken water depletion into account, 
these planets contain about twenty times as much water 
as those formed in sim i lar out er giant planet simulations 
(|Ravmond et al.ll2004 l2006d) . These planets are likely 
to be " water worlds", cov ered in kilometers-deep global 
oceans (jLeger et al.ll2004 ). 

The variations in the evolutionary process for each of 
the three models examined here demonstrate the impor- 
tance of both the configuration of giant planets in the sys- 
tem as well as the presence and duration of the gaseous 
disk in dictating the final parameters of the terrestrial 



planets formed. Simulations run without including the 
gaseous disk (JSN) showed very little planetesimal ac- 
cretion during migration and almost no collisions after 
migration due to the strong dynamical excitation of the 
disk (see Figure[7|). Though embryo eccentricities and in- 
clinations decreased due to dynamical friction, the mean 
embryo composition remained similar to the initial value 
and final terrestrial planet masses were in the lower end 
of the mass range seen in the Solar System (see Figure 
lll[) . In contrast, embryos in simulations including the 
presence of a dissipating gaseous disk (JD and JSD) ac- 
creted planetesimals and other embryos from both the 
inner system ( during Stage 1) and the outer system 
(during Stage 2). This resulted in both higher planet 
masses and water mass fractions, though this effect was 
slightly mitigated in the JSD simulations due to more 
rapid clearing in the outer system. Eccentricities and in- 
clinations are also affected by the presence of an outer 
giant planet; in the JD simulations dynamical excitation 
continues unabated during Stage 3, while in the JSD sim- 
ulations excited bodies are scattered by the Saturn-mass 
planet and the remaining bodies therefore remain dy- 
namically relaxed. 

5.2. Implications for Our Own System 

If planetary systems which suffer the migration of a gi- 
ant planet to small distances can eventually form terres- 
trial planets similar to those in our own system, and the 
migration of young giant planets is a common result of 
interactions with the gaseous disk, then it is appropriate 
to consider the possibility that our own planetary sys- 
tem could have formed earlier generations of giant plan- 
ets prior to those in the outer Solar System. The systems 
simulated here undergo extensive radial mixing of large 
bodies throughout the inner and outer system through 
scattering, which could leave a elemental signature on 
the final population of both planetary and sub-planetary 
bodies. However, radial mixing of dust and smaller bod- 
ies is also thought to occur through a variety of other 
processes in the pla net formation r egion besides plane- 
tary scattering (see IWooden et al.l ((2005) for a detailed 
discussion), and it is clear from the final abundances of 
planets in our simulations and others that any signature 
of the scattering of massive bodies early in the formation 
history will be quickly erased through accretion in the 
inner system due to continued mixing. This is less true 
in the outer system, where long dynamical timescales 
and low solid and gas densities make accretion and ra- 
dial transport less likely. Therefore one of the primary 
constraints on giant planet formation and migration the- 
ories is the scattered population of Kuiper Belt objects 
(KBO's) and Oort Cloud comets - both their orbital and 
compositional characteristics. 

Though conclusive detection of anomalous cometary 
and asteroidal abundances is naturally difficult and am- 
biguous, the detection of comet-like bodies with orbital 
characteristics suggesting an origin in the Oort Cloud but 
with physical appearances similar to asteroids could be 
considered a strong indicator of the scattering of mas- 
sive inner-disk bodies by a migrating giant planet. In 
fact, such bodies (known as Damocloids) have recently 
been discovered (jAsher et al.lll994 lJewittll2005ft . but it 
is unclear whether any of these bodies are truly aster- 
oidal in composition, and even if they are it has been 
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shown that the stochastic nature of late-stage scatter- 
ing by Jupiter could also produce s mall numbers of Port 
Clou d comets of asteroidal nature (jWeissman &; Levisonl 
Bug) , a more conclusive sign of giant planet migration 
would be a classical KBO with a composition primar- 
ily composed of refractory materials; this would imply 
the re-circularization of a scattered inner-disk object, 
which would most likely only be possible in the pres- 
ence of damping by significant amounts of gas or dust 
for long timescales. In-depth analysis of the meteoritic 
and cometary record and other signatures of long-term 
planetary dynamics for evidence of early giant planet mi- 
gration is clearly beyond the scope of this paper, and 
would most likely be inconclusive considering the many 
factors involved. We must wait for the detection and 
accurate analysis of many more KBO's and Oort Cloud 
comets before we can place meaningful constraints on the 
possibility of Type II giant planet migration in our own 
system. 

5.3. Future Work 

In this study we have endeavored to utilize the most 
concrete data and plausible scenarios to develop our ini- 
tial conditions for each model, while limiting the number 
of free parameters and ill-constrained variables. How- 
ever, there is naturally a range of parameter space in the 
formation of terrestrial planets in systems with a migrat- 
ing giant planet that we could not explore. 

We have only considered systems with one migrating 
giant planet. The observed distribution of giant planet 
orbital parameters can be roughly reproduced by the 
comb ination of migratio n and gravitational scattering 
(e.g.. lAdams fe Laughlinl (|2003l ). Thus, in certain cases 
there may be multiple giant planet s mi grating toge ther. 
Simulations bv lBrvden et all (|2000f ) and lKlevi i|2000f ) sug- 
gest that if the planets attain their full mass while on 
closely-spaced orbits, the intervening gas between them 
would cause convergent migration, leading to resonance 
capture or scattering. However, if one planet forms suf- 
ficiently early to migrate inwards, a second planet could 
form beyond the evacuation zone and migrate as well. 
In this case gravitational scattering would be greatly en- 
hanced because inner disk material would interact with 
both giant planets and would probably be scattered to 
larger orbital distances (as occurs for some bodies in 
model JSD). However, in most cases the dynamic insta- 
bilities of the two massive planets would likely out-weigh 
the evolution of the smaller bodies, leading to chaotic 
scatteri ng and removal of small bodies similar as out- 
lined in lVeras fc Armita ge ( 2005). If a combined migra- 
tion leaves one planet in a small orbit within 0.5 AU 



and a second planet beyond 2 AU, on e would expect a 
hybrid scenario between this work and iRavmondl (|2006l ) 
with terrestrial planets forming in the stable region be- 
tween the two giant planets, but more work is needed to 
confirm this. 

Additional uncertainty surrounds the role of the 
gaseous disk. In this work we bracket the potential range 
of disk masses with the JSD and JSN models. As stated, 
the differences in the damping of scattering energies and 
the orbital decay of planetesimal orbits greatly impacts 
the final characteristics of the terrestrial planets that are 
formed. To fully investigate the full range of gas disk 
properties and planetesimal size ranges would require an 
unreasonably large number of simulation runs. Our sim- 
ulations were limited by computational resources. We 
chose fixed values for several parameters such as the mi- 
gration rate and disk properties, and chose to integrate 
these simulations for a long time rather t han explore the 
effects of different parameters. (Note that lFogg fe Nelson] 
(2005) explored the effects of the disk's mass distribu- 
tion.) In addition, we did not include the effects of 
collisional fragmentation, which are certainly impo rtant 
for both planetesimal- f e.g.JBenz fe Asphaud (|1999f ) and 
embryo-scale collisions (|Asphaug et al.ll2006t )). More ef- 
fects relating to the gas disk also remain to be modele d, 
including tidal gas drag (e.g.. iKominami fe Ida! ((2002)), 
which acts preferentially on larger bodies. Clearly, more 
work remains to include all of the relevant physics. De- 
spite these limitations, our simulations are among the 
first to realistically address this problem and examine the 
final results of terrestrial planet formation under these 
condit i ons (also see earlier w o rk: iMandell fc Sigurdsson! 
(|2T)0l):lFbggJ L Nelsonl ((2005h : iZhou et all (|2005h : Paper 
1; iFogg fc Nelson! (|2006f )). Hopefully this work will en- 
courage future studies to expand and improve the models 
explored here, both in our understanding of the relevant 
physical conditions and our ability to realistically model 
the myriad and complex forces at work in the formation 
of terrestrial planetary systems. 
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Fig. 1. — Initia l value s for the mass and water mass fraction of protoplanetary disk used for each simulation. The solid disk is modeled 
after the Hayashi (1981) MMSN, augmented by a factor of 2.2. Embryos are spaced randomly by A=5-10 mutual Hill radii, and embryo 
masses increase with radial distance r as M em j, oc A 3 / 2 r 3 / 4 . Water and iron contents are based on measurements from Solar System bodies 
l|Abe et al.ll2000l : ILodders fc Feglevl[T998T l. 
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Fig. 2. — Evolution of a sample simulation run using the JD model (a single migrating Jupiter-mass planet with gas drag included). 
Radial distance is plotted on the x-axis, with eccentricity on the y-axis. Relative inclination is indicated through error-bars on each point, 
and water content is indicated by the color (a reference bar is located below the plot). The size of each point is proportional to M 1 ' 3 , with 
the size of the inner black region representing the iron mass fraction. The system forms a 'hot Earth' and a massiv e terrestrial plan e t in the 
Habitable Zone. Water contents are high through radial mixing and migration of planetesimals. Reprinted from|Ravmond ct al. (2006b) 
(Paper I). 
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Fig. 3. — Evolution of a sample simulation run using the JSD model (a migrating Jupiter-mass planet with a second Saturn-mass outer 
planet and gas drag included). The evolution on the inner system is similar to the JD model, but the outer system is cleared rapidly and 
less water-rich material flows into the inner system. 
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Fig. 4. — Evolution of a sample simulation run using the JSN model (a migrating Jupiter-mass planet with a second Saturn-mass outer 
planet and no gas drag). The protoplanetary disk is highly excited, with most the the mass being lost to ejection. Over long time scales 
dynamical friction cools the disk, and low-mass planets with lower water contents form. 
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Fig. 5. — The average total mass of the three models (in Mj) versus time. The total mass is divided into three zones: the inner disk 
(a < 5.2AU), the intermediate disk (5.2AU < a < 9.5AU), and the outer disk (a > 9.5AU). Additionally, a dashed line is placed at the 
embryo / planetesimal mass boundary - below the line is the mass from embryos, above it is the mass from planctcsimals. Differences can 
be seen between the three models, specifically in the loss of mass without the presence of gas and the rapid evolution of the outer disk with 
the presence of a second giant planet. 
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Fig. 6. — Illustration of the loss processes for embryos and planetesimals over time for the three models. Accretion happens primarily 
during migration, with additional accretion occurring as the gas disappears at 10 Myr. Simulations with gas drag show almost complete 
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Fig. 7. — Number of bodies accreted as a function of time, with the radial location of the impact denoted by color, for the three different 
models. Accretion occurs almost solely within 1.5 AU during migration, with additional accretion occurring at a range of distances as the 
gas disk dissipates. Accretion is almost non-existent when the gas is dense enough to damp down orbital excitation. 
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are shepherded in resonances to become 'hot Earths'. Other embryos are scattered outward to become outer terrestrial planets. Multiple 
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Fig. 9. — Plot of mass versus time (top) and water mass fraction versus time (bottom) for three indicative bodies in a sample JD 
simulation. The body that settles at 0.12 AU (a 'hot Earth') accretes very rapidly near the end of migration, while the outer terrestrial 
planets accrete material gradually over longer timescales. For the outer terrestrial planets, most of their water is accreted at two distinct 
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Fig. 10. — Final configurations for 'hot Earths' in all simulations with gas drag (JD and JSD) before instability sets in when the gas disk 
dissipates. Resonance positions for a 'hot Jupiter' at 0.25 AU are marked. The size of each body corresponds to its mass, with the orbital 
range (i.e. eccentricity) indicated below. It is clear that most 'hot Earths' reside just inside resonances, due to the combined shepherding 
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Fig. 11. — Final configurations for all the simulations after 200 Myr of integration. The size of each body indicates mass and the color 
indicates water content (with the Solar System plotted on the bottom for reference). The orbital range (i.e. eccentricity) for each body 
is indicated by the horizontal error-bar. It is clear that systems with long-lived gas disks will form more massive, water-rich planets in or 
near the Habitable Zone, but that planets will form even in systems without gas. Additionally, dense gas in the inner system will lead to 
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Fig. 12. — The known extrasolar planets that are likely to have formed terrest rial planet s in th e Habitable Zone. The solid line shows 
the outer giant planet semi-major axis limit for habitable planet formation from Raymond (2006), and the dashed line is the inner giant 
planet semi-major axis limit found from the current simulations. Filled circles indicate the known giant planets that satisfy our criteria 
(eccentricity less than 0.1 and semi-major axis within the derived limits) and may therefore harbor a terrestrial planet in the Habitable 
Zone; open circles are unlikely candidates for habitab ility. Note that not all p lanets above or below our limits qualify, due to the eccentricity 
limit. The Habitable Zone is shaded. Modified from Raymond ct al. (2006b) (Paper I). 
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TABLE 1 

Embryo/Planetesimal Properties at the End of Stage 1 





n 1 

Zone 1 


JD 


JSD 


JSN 




Z.l 


7.60 M / 0.31 M e 


7.97 M ffl / 0.53M e 


1.90 M ffi / 1.43 M e 


Tot Mass- Em /PI 


Z.2 


6.18 M® / 1.42 Mj) 


3 24Mib / 70 Mm 


4 18 Mm / 1 37 Mm 




Z.3 


2.77 Mm 1 0.07 Mm 


4.93 Mm / 0.50 Mm 


3.36 M e / 1.79 Mm 


Em. M (Mm) 2 


Z.l 


0.45 (0.07 - 4.35) 


0.46 (0.08 - 4.90) 


0.18 (0.05 - 0.54) 


Z.2 


0.42 (0.08 - 1.23) 


0.29 (0.09 - 1.10) 


0.39 (0.05 - 1.05) 




Z.3 


0.28 (0.09 - 0.88) 


0.39 (0.10 - 1.10) 


0.32 (0.03 - 1.23) 




Z.l 


-2.09 (-5.00 - -1.26) 


-1.74 (-5.00 - -0.30) 


-1.75 (-5.00 - -0.30) 


Em. log(W.M.F.) 2 


Z.2 


-0.56 (-5.00 - -0.30) 


-0.72 (-5.00 - -0.30) 


-0.59 (-5.00 - -0.30) 




Z.3 


-0.85 (-5.00 - -0.30) 


-0.60 (-5.00 - -0.30) 


-0.75 (-5.00 - -0.30) 




Z.l 


-0.52 (-1.13 - -0.37) 


-0.54 (-1.18 - -0.40) 


-0.51 (-1.01 - -0.34) 


Em. log(Fe M.F.) 2 


Z.2 


-0 76 (A 23 - -0 34") 


-0 74 (-1 23 - -0 37*1 


-0 74 (-1 23 - -0 38") 




Z.3 


-0.80 (-1.23 - -0.49) 


-0.81 (-1.23 - -0.39) 


-0.70 (-1.23 - -0.37) 




Z.l 


0.44 (0.01 - 0.84) 


0.48 (0.01 - 0.87) 


0.61 (0.02 - 0.91) 


Em. Eccen. 2 


Z.2 


37 (0 01 - 92") 


42 (0 06 - 86") 

KJ ."-t£s \KJ.\JKJ VJ > (J\J J 


44 (0 04 - 94 ) 




Z.3 


0.75 (0.41 - 0.97) 


0.62 (0.09 - 0.99) 


0.71 (0.06 - 0.98) 




Z.l 


6.62 (0.22 - 36.65) 


8.49 (0.20 - 29.28) 


26.34 (0.51 - 80.41) 


Em. Inclin.(°) 2 


Z.2 


5.25 (0.44 - 28.00) 


12.89 (0.90 - 51.13) 


19.79 (1.26 - 70.18) 


Z.3 


4.69 (0.64 - 13.56) 


11.79 (1.23 - 34.50) 


17.09 (1.40 - 46.70) 




Z.l 


0.00 (0.00 - 0.10) 


0.01 (0.00 - 0.19) 


0.54 (0.06 - 0.94) 


PI. Ecccn. 2 


Z.2 


0.09 (0.00 - 0.23) 


0.08 (0.01 - 0.40) 


0.41 (0.01 - 0.97) 




Z.3 


0.11 (0.01 - 0.21) 


0.31 (0.02 - 0.79) 


0.60 (0.02 - 0.99) 




Z.l 


0.08 (0.01 - 1.86) 


0.50 (0.00 - 5.51) 


28.65 (0.19 - 84.56) 


PI. Inclin.(°) 2 


Z.2 


2.58 (0.02 - 10.52) 


3.67 (0.23 - 18.00) 


17.93 (0.33 - 80.30) 




Z.3 


2.36 (0.61 - 5.00) 


5.69 (0.35 - 32.37) 


14.94 (0.30 - 81.01) 



Radial zones arc delineated as follows: Z.l (a<5.2 AU), Z.2 (5.2<a<9.5), and Z.3 (9.5<a) 2 Values represent the mean of the bodies in each zone, 
with the range of values in parentheses. 
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TABLE 2 

Embryo/Planetesimal Properties at the End of Stage 2 





Zone 1 


JD 


JSD 


JSN 


Tot. Mass: Em/Pi 


Z.l 
Z.2 
Z.3 


9.25 M© / 1.01 M e 
5.12 M© / 0.02 M® 
2.65M© / 0.00M© 


5.77 M® / 0.51 M® 
0.63 M© / 0.01 M® 
0.26 M® / 0.44 M® 


1.34 M® / 0.82 M® 
0.67M® / 0.19 M© 
1.04 M® / 0.24 M® 


Em. M (M e ) 2 


Z.l 
Z.2 
Z.3 


0.55 (0.09 - 4.59) 
0.46 (0.12 - 1.59) 
0.32 (0.11 - 0.93) 


0.40 (0.08 - 4.44) 
0.42 (0.14 - 0.95) 
0.52 (0.15 - 0.90) 


0.24 (0.11 - 0.81) 
0.38 (0.14 - 0.60) 
0.52 (0.12 - 1.05) 


Em. log(W.M.F.) 2 


Z.l 
Z.2 
Z.3 


-1.03 (-5.00 - -0.30) 
-0.54 (-4.41 - -0.30) 
-0.72 (-5.00 - -0.30) 


-1.03 (-5.00 - -0.30) 
-0.57 (-2.74 - -0.30) 
-0.56 (-1.30 - -0.30) 


-1.26 (-5.00 - -0.30) 
-0.53 (-5.00 - -0.30) 
-0.50 (-5.00 - -0.30) 


Em. log(Fc M.F.) 2 


Z.l 
Z.2 
Z.3 


-0.59 (-1.13 - -0.42) 
-0.81 (-1.23 - -0.48) 
-0.86 (-1.23 - -0.49) 


-0.61 (-1.20 - -0.40) 
-0.83 (-1.13 - -0.55) 
-0.99 (-1.18 - -0.86) 


-0.56 (-1.01 - -0.41) 
-0.75 (-0.96 - -0.48) 
-0.79 (-1.22 - -0.49) 


Em. Eccen. 2 


Z.l 
Z.2 
Z.3 


0.13 (0.00 - 0.55) 
0.29 (0.02 - 0.86) 
0.53 (0.17 - 0.93) 


0.07 (0.00 - 0.64) 
0.17 (0.01 - 0.40) 
0.73 (0.73 - 0.74) 


0.34 (0.01 - 0.79) 
0.40 (0.11 - 0.85) 
0.42 (0.15 - 0.80) 


Em. Inclin.(°) 2 


Z.l 
Z.2 
Z.3 


2.63 (0.03 - 16.09) 
6.85 (0.72 - 21.31) 
6.77 (0.83 - 16.80) 


2.07 (0.05 - 26.80) 
5.53 (0.91 - 13.54) 
21.40 (15.11 - 30.44) 


20.87 (3.00 - 48.33) 
32.45 (7.42 - 55.09) 
22.97 (5.55 - 60.20) 


PI. Ecccn. 2 


Z.l 
Z.2 
Z.3 


0.06 (0.00 - 0.31) 
0.16 (0.05 - 0.25) 
0.21 (0.18 - 0.23) 


0.03 (0.00 - 0.31) 
0.05 (0.02 - 0.11) 
0.05 (0.00 - 0.29) 


0.47 (0.04 - 0.85) 
0.41 (0.03 - 0.92) 
0.57 (0.06 - 0.97) 


PI. Inclin.(°) 2 


Z.l 
Z.2 
Z.3 


1.48 (0.01 - 12.63) 
8.03 (1.69 - 11.38) 
6.54 (4.21 - 9.10) 


1.18 (0.05 - 6.29) 
1.73 (1.41 - 2.56) 
1.82 (0.24 - 18.63) 


27.68 (0.74 - 69.72) 
27.78 (3.14 - 66.70) 
27.09 (1.04 - 71.22) 


1 Radial zones are delineated as follows: Z.l (a<5.2 AU), Z.2 (5.2<a<9.5), 
with the range of values in parentheses. 


and Z.3 (9.5<a) 2 Values represent the mean 


of the bodies in each zone 
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TABLE 3 

Final Embryo/Planetesimal Properties 





Zone 1 


JD 


JSD 


JSN 


Tot. Mass: Em/Pi 


Z.l 
Z.2 
Z.3 


5.52 M ffi / 0.05 M e 
2.71 M ffi / 0.00 M ffi 
4.83 M ffi / 0.02 M ffi 


4.08 M e / 0.05 Me 
0.26 Me / 0.00 M ffi 
0.00 M ffi / 0.35 M e 


1.17M ffi / 0.19 Me 
0.00 M ffi / 0.01 M e 
0.51 M ffi / 0.04M e 


Em. M (M e ) 2 


Z.l 
Z.2 
Z.3 


1.03 (0.12 - 4.15) 
0.48 (0.14 - 1.59) 
0.32 (0.09 - 1.24) 


0.58 (0.08 - 1.77) 
0.35 (0.17 - 0.46) 
0.00 (0.00 - 0.00) 


0.29 (0.11 - 0.82) 
0.00 (0.00 - 0.00) 
1.03 (1.01 - 1.05) 


Em. log(W.M.F.) 2 


Z.l 
Z.2 
Z.3 


-0.68 (-2.93 - -0.30) 
-0.56 (-5.00 - -0.30) 
-0.67 (-5.00 - -0.30) 


-0.87 (-5.00 - -0.30) 
-0.71 (-1.46 - -0.30) 
0.00 (0.00 - 0.00) 


-1.01 (-5.00 - -0.30) 

0.00 (0.00 - 0.00) 
-0.30 (-0.30 - -0.30) 


Em. log(Fc M.F.) 2 


Z.l 
Z.2 
Z.3 


-0.65 (-0.96 - -0.50) 
-0.81 (-1.21 - -0.54) 
-0.81 (-1.23 - -0.48) 


-0.65 (-1.19 - -0.44) 
-0.69 (-0.88 - -0.59) 
0.00 (0.00 - 0.00) 


-0.55 (-0.94 - -0.42) 

0.00 (0.00 - 0.00) 
-0.91 (-0.93 - -0.90) 


Em. Eccen. 2 


Z.l 
Z.2 
Z.3 


0.29 (0.01 - 0.74) 
0.30 (0.08 - 0.58) 
0.39 (0.06 - 0.95) 


0.17 (0.00 - 0.57) 
0.15 (0.12 - 0.35) 
0.00 (0.00 - 0.00) 


0.24 (0.03 - 0.42) 
0.00 (0.00 - 0.00) 
0.19 (0.13 - 0.24) 


Em. Inclin.(°) 2 


Z.l 
Z.2 
Z.3 


14.40 (0.69 - 32.24) 
18.46 (4.73 - 42.86) 
19.76 (2.71 - 50.25) 


6.66 (0.35 - 21.80) 
13.66 (6.50 - 19.22) 
0.00 (0.00 - 0.00) 


17.27 (3.71 - 28.34) 
0.00 (0.00 - 0.00) 
3.21 (2.74 - 3.82) 


PI. Ecccn. 2 


Z.l 
Z.2 
Z.3 


0.13 (0.05 - 0.34) 
0.67 (0.67 - 0.67) 
0.69 (0.40 - 0.90) 


0.44 (0.03 - 0.89) 
0.37 (0.37 - 0.37) 
0.08 (0.01 - 0.26) 


0.40 (0.06 - 0.81) 
0.43 (0.17 - 0.73) 
0.28 (0.06 - 0.60) 


PI. Inclin.(°) 2 


Z.l 
Z.2 
Z.3 


8.80 (2.28 - 50.79) 
93.39 (93.39 - 93.39) 
33.94 (12.39 - 92.93) 


31.25 (2.02 - 61.23) 
43.84 (43.84 - 43.84) 
4.07 (0.21 - 27.00) 


31.65 (2.02 - 59.38) 
31.21 (8.42 - 48.70) 
18.87 (4.03 - 42.56) 



Radial zones arc delineated as follows: Z.l (a<5.2 AU), Z.2 (5.2<a<9.5), and Z.3 (9.5<a) 2 Values represent the mean of the bodies in each zone, 
with the range of values in parentheses. 
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TABLE 4 

Properties of planets formed in JD Simulations 1 



Simulation 


a(AU) 


e 2 


in 


M(M ffi ) 


W.M.F. 


Fe M.F. 


JD-1 


0.10 


0.01 


0.6 


323.42 


- 


- 




0.64 


0.33 


1.5 


1.14 


7.72 x 10~ 2 


0.27 




1.80 


0.33 


5.5 


2.14 


1.13 x lO" 1 


0.27 




3.95 


0.47 


7.2 


0.90 


5.00 x lO" 1 


0.11 




4.01 


0.46 


24.4 


0.24 


1.00 x 10" 5 


0.29 


JD-3 


0.14 


0.01 


0.5 


326.45 




— 




1.74 


0.46 


38.0 


0.32 


5.78 x 10~ 2 


0.29 




2.27 


0.48 


22.9 


0.36 


1.01 x 10 _1 


0.27 




3.80 


0.27 


9.0 


1.25 


4.05 x 10 _1 


0.16 


JD-4 


0.14 


0.02 


0.3 


324.17 










0.31 


0.06 


1.4 


1.63 


4.60 x 10~ 2 


0.31 




1.19 


0.12 


8.0 


0.92 


7.64 x 10~ 2 


0.27 




2.47 


0.59 


42.7 


0.30 


1.52 x 10 _1 


0.24 




2.86 


0.21 


10.5 


1.06 


7.08 x 10~ 2 


0.20 




3.33 


0.31 


24.0 


0.23 


9.23 x 10~ 2 


0.28 




3.40 


0.19 


12.5 


0.83 


5.00 x 10 _1 


0.11 




3.67 


0.21 


16.3 


0.14 


7.39 x 10~ 2 


0.26 




3.80 


0.69 


45.7 


0.14 


2.28 x 10 _1 


0.22 




3.92 


0.23 


11.2 


0.39 


9.03 x 10~ 2 


0.26 




4.60 


0.32 


10.3 


0.46 


1.79 x lO" 1 


0.17 


JD-5 


0.12 


0.01 


0.7 


4.15 


2.34 x 10~ 2 


0.25 




0.21 


0.01 


0.1 


320.92 








0.91 


0.08 


5.0 


3.05 


9.23 x 10~ 2 


0.27 




2.98 


0.33 


13.3 


0.55 


1.36 x lO" 1 


0.26 




3.38 


0.22 


11.9 


0.32 


1.15 x 10" 3 


0.31 




4.32 


0.51 


27.7 


0.12 


3.77 x 10" 3 


0.27 




4.78 


0.27 


14.9 


0.84 


5.00 x 10" 1 


0.13 




4.68 


0.14 


13.5 


0.72 


5.00 x 10- 1 


0.11 



Planets are defined to be > 0.1 and inside 5 AU. Close-in giant planets are shown in bold. Orbital elements are averaged over the last Myr 
of the simulation. 
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TABLE 5 

Properties of planets formed in JSD Simulations 



Simulation 


a(AU) 


e 




M(M®) 


W.M.F. 


Fe M.F. 


JSD-1 


0.08 


0.00 


0.2 


325.57 




— 




0.19 


0.01 


0.5 


0.46 


1.12 x 10 _1 


0.25 




0.86 


0.12 


4.1 


0.46 


4.76 x 10~ 2 


0.20 




1.62 


0.07 


2.0 


1.46 


4.34 x 10~ 2 


0.23 




2.32 


0.22 


3.4 


0.46 


2.59 x 10~ 2 


0.30 




3.38 


0.08 


2.6 


0.29 


7.68 x 10~ 2 


0.23 




4.57 


0.12 


4.6 


0.46 


5.00 x lO" 1 


0.13 


JSD-2 


0.06 


0.00 


2.3 


325.67 




— 




0.14 


0.02 


2.3 


0.38 


9.99 x 10 -2 


0.13 




0.26 


0.00 


2.3 


0.14 


5.55 x 10~ 2 


0.28 




0.75 


0.29 


17.2 


0.41 


1.43 x 10 _1 


0.23 




1.42 


0.18 


7.1 


1.77 


1.24 x 10 _1 


0.27 




3.09 


0.23 


5.1 


0.95 


5.00 x 10 _1 


0.12 




4.06 


0.32 


8.9 


0.60 


5.59 x 10~ 2 


0.23 


JSD-3 


0.07 


0.01 


0.7 


1.24 


3.25 X 10 


0.25 




0.11 


0.00 


0.6 


3.23 


1.92 x 10~ 2 


0.31 




0.19 


0.00 


0.6 


320.21 


- 


- 




0.41 


0.03 


1.4 


0.12 


1.00 x 10~ 5 


0.35 




0.51 


0.01 


1.0 


0.68 


3.67 x 10~ 2 


0.32 




0.61 


0.01 


1.1 


0.74 


4.14 x 10~ 2 


0.32 




1.17 


0.02 


0.8 


0.24 


7.57 x 10~ 2 


0.23 




2.34 


0.12 


4.1 


0.19 


3.70 x 10~ 2 


0.28 




2.63 


0.17 


6.5 


0.13 


2.90 x 10~ 2 


0.31 




3.67 


0.10 


4.8 


0.33 


6.68 x 10~ 2 


0.16 




3.93 


0.18 


15.0 


0.17 


6.81 x 10~ 2 


0.20 




4.49 


0.17 


2.7 


0.61 


5.00 x lO" 1 


0.12 


JSD-4 


0.05 


0.00 


2.1 


324.82 








0.25 


0.01 


2.1 


1.11 


5.20 x 10~ 2 


0.20 




0.55 


0.14 


8.1 


0.77 


1.05 x lO" 1 


0.19 




0.92 


0.14 


6.8 


0.36 


1.37 x 10- 1 


0.27 




1.26 


0.07 


7.1 


0.63 


1.16 x lO" 1 


0.21 




2.65 


0.20 


19.8 


0.11 


3.33 x 10" 3 


0.27 




3.02 


0.04 


2.0 


1.32 


3.25 x 10- 1 


0.17 
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TABLE 6 

Properties of planets formed in JSN Simulations 



Simulation 


a(AU) 


e 




M(M®) 


W.M.F. 


Fe M.F. 


JSN-2 


0.24 


0.07 


1.1 


318.85 








1.13 


0.35 


24.9 


0.19 


1.84xl0~ 3 


0.30 




1.32 


0.20 


17.0 


0.28 


1.26X10 -3 


0.35 




3.98 


0.06 


10.2 


0.53 


4.93X10- 1 


0.13 


JSN-3 


0.05 


0.37 


29.3 


0.22 


4.99X10- 1 


0.11 




0.24 


0.05 


3.7 


318.82 








1.70 


0.25 


3.8 


0.24 


4.85xl0~ 2 


0.24 




2.26 


0.04 


7.0 


0.18 


lxlO- 5 


0.30 




4.04 


0.40 


21.2 


0.22 


lxl0~ 5 


0.29 




4.77 


0.18 


14.0 


0.54 


5.00X10- 1 


0.14 




22.70 


0.25 


3.8 


1.01 


5.00X10- 1 


0.12 


JSN-4 


0.24 


0.04 


1.7 


319.29 








0.96 


0.07 


17.3 


0.22 


2.46X10 -3 


0.28 




2.76 


0.04 


28.9 


0.30 


2.13xl0~ 3 


0.26 




4.50 


0.22 


12.4 


0.82 


5.00X10- 1 


0.11 




18.53 


0.13 


3.6 


1.05 


5.00X10- 1 


0.13 


JSN-5 


0.24 


0.04 


2.0 


319.10 








1.10 


0.21 


12.0 


0.40 


lxlO" 5 


0.36 




1.15 


0.31 


24.6 


0.13 


5.44xl0" 3 


0.29 




2.68 


0.28 


13.3 


0.17 


2.02X10" 3 


0.38 




3.39 


0.47 


23.4 


0.11 


lxlO" 5 


0.31 




3.59 


0.24 


20.0 


0.12 


lxlO" 5 


0.32 
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TABLE 7 

Giant Planet Semi-Major Axis Limits for Potentially Habitable Systems 



M* (M n ) 


Sp. Type 1 


Hab Zone (AU) 2 


Inner Limit (AU) 


Outer Limit (AU) 


0.1 


M6 


0.024 - 0.045 


0.015 


0.075 


0.4 


M3 


0.10 - 0.19 


0.06 


0.32 


0.7 


K6 


0.28 - 0.52 


0.17 


0.87 


1.0 


G2 


0.8 - 1.5 


0.5 


2.5 


1.3 


F8 


2.3 - 4.3 


1.45 


7.2 


1.6 


F0 


6.5 - 12.3 


4.1 


20.5 


2.0 


A5 


25 - 47 


15.7 


78.3 


1 Spectral types 
age-dependent . ^ 


from Table 8.1 oflReid & Hawlevl J2000I) and Appendix E from ICarroll & Ostliel 119961). Spectral types 
Habitable Zones scaled by L*/ 2 to 0.8 - 1.5 AU for a solar-mass star. 


of low-mass stars are 
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TABLE 8 

Potentially Habitable Exoplanet Systems 



System 


M* (M ) 


[Fc/H] 


M p ; (Mj) 


a (AU) 


e 


HZ 


Om\F 0^ 071T, 

KJ \_T LjUj U'J U I J_J_j 


0.13 




0.9 


1.800 




u.uo u.uu 




0.22 




0.02 


2.100 




n or n i n 

U . UU U. J-U 


VjJ OO -L 


0.31 


-0.25 


0.052 


0.041 


0.00 


n ns n 1 4 

u .uo U. ±1 


VJT 1_J i_J At 'J O IV J. VAUO 


0.36 




1.5 


3.000 




fl OQ-0 1 7 

U.Uf U. J- I 


UF) 41 004 A 


0.40 


0.16 


17.892 


0.018 


0.08 


n 1 n n i q 

u . _Lu u. J. y 


nru,F o^ i^qt, 


0.49 




0.04 


2.800 




14 9^ 

U . -L rr U . ^ U 


un ^007^ 


0.70 


0.08 


0.623 


0.039 


0.00 


n 98 n ^9 

U.ZrO U. OA 


UF) 978Q4 

L1LJ Zj 1 Of rt 


0.75 


0.3 


0.618 


0.122 


0.05 


u.oo u. u.^ 


HD 1 1 4386 


0.76 


0.004 


1.343 


1.714 


0.23 


^4-0 64 


V\T) 1 ^44^ 


0.77 


-0.27 


3.9 


0.113 


0.04 


U.OU U. UU 


OdTF TR 1 1 3 


0.77 


0.14 


1.35 


0.023 


0.00 


n ^ n fifi 

U.OU U. UU 


HD 1 1 1 232 


0.78 


-0.36 


6.803 


1.975 


0.20 


fl ^7-0 fiQ 




0.80 


0.11 


0.385 


0.036 


0.00 


n n 74 

u.oy u. / 1 


HD 1Q22fi^ 

11J_/ 1l/^ZiUU 


0.81 


0.05 


0.641 


0.153 


0.05 


0.41-0.76 


ocit.f tr 1 1 1 


0.81 


0.12 


0.52 


0.047 


0.00 


n 41 n 7fi 

U.I J. U. 1 u 




0.82 


-0.03 


1.058 


3.377 


0.25 


49 7Q 
u.i^ u. i y 


HD 18Q7^S 

-1 -I J—* _LOy I UO 


0.82 


-0.03 


1.152 


0.031 


0.00 


fl 42-0 7Q 

U.rlii u. i y 


T4D 1 ^0^99 


0.88 


0.006 


1.088 


0.091 


0.02 


^9 


TrES-1 


0.89 


0.001 


0.759 


0.039 


0.00 


0.54-1.01 


T4D 4^08 


0.90 


-0.31 


0.047 


0.118 


0.00 


^fi 1 0^1 

U.UU l.UU 


OO V_;I1U c 


0.91 


0.31 


0.038 


0.038 


0.09 


^8 1 HQ 


^ Pnr h 

UU Vyllt, u 






0.833 


0.114 


0.01 


1 OQ 

u .oo J. . uy 


55 Cue c 






n K7 

U. lO I 


n 9^8 

U. ZrOO 


n 07 


a w i no, 
u.oo-j-.uy 


UU V^Ut., Ll 






3.887 


5.964 


0.09 


^8 1 OQ 
u. uo J. . uy 


T4D 4fi^7^ 
njj iuo / u 


0.92 


0.24 


0.226 


0.040 


0.06 


fiO 1 1 ^ 

U.UU -L . -LO 


T4D 9fi^8 


0.93 


0.16 


0.477 


0.044 


0.00 


0.62-1.17 


HD 1 021 95 


0.93 


-0.09 


0.492 


0.049 


0.06 


0.62-1.17 


UP) 1fizLQ99 


0.94 


0.17 


0.36 


2.110 


0.05 


fi^i 1 91 

U.UU -L . -L 


T4D 8Q307 


1.00 


-0.16 


2.601 


3.945 


0.01 


80 1 ^0 

U . OU -L . UU 


HD 8^443 


1.00 


0.36 


0.398 


0.041 


0.01 


80-1 ^0 

U . O U l.JU 


1_\,1±U V^Ill., J_> 


1.00 


-0.20 


1.092 


0.229 


0.06 


80 1 ^0 

U . OU -L . UU 


rd -1 n 3i fifi 

III/ 1 \ I OIUU 


1.01 


0.38 


0.458 


0.045 


0.02 


8*3-1 ^ 

U.OU J. . UU 


UP) 70fi49 


1.05 


0.16 


1.97 


3.230 


0.03 


Qfi 1 7Q 
u. yu J. . i y 


up) 7^9^fi 

L1LJ 1 0.£UU 


1.05 


0.26 


1.867 


0.037 


0.03 


Qfi 1 7Q 
u. yu i. (y 


T4D 91 9^01 


1.05 


-0.18 


0.396 


0.034 


0.00 


Qfi 1 7Q 
u. yu J. . / y 


T4D 4Qfi74 

J.1.L/ Try U 1 rfc 


1.06 


0.31 


0.105 


0.058 


0.09 


QQ 1 8fi 
u.yy _L .ou 


fm 195019 

11J_/ _L£7UU-Ly 


1.07 


0.07 


3.681 


0.139 


0.01 


1.03-1.92 


T4D 1 871 9^ 

LlLy _LO I _L^O 


1.08 


0.12 


0.527 


0.043 


0.02 


1 Ofi 1 QQ 
J. . uu J. . y y 


^1 Ppcr 

01 r eg 


1.09 


0.20 


0.472 


0.053 


0.01 


1.10-2.07 


T4D 107148 


1.12 


0.31 


0.21 


0.269 


0.05 


1 93 9 ^0 

J_ .ZiU Zi.OU 


un 7fi700 


1.13 


0.35 


0.232 


0.051 


0.09 


1 97 9 ^8 


HD 20Q4 r SS 


1.14 


0.014 


0.689 


0.047 


0.00 


1.31-2.47 


OOT,F TR ID 


1.17 


0.12 


0.63 


0.042 


0.00 


1.46-2.74 


OnT,F TR 5fi 


1.17 




1.24 


0.023 




1.46-2.74 


HD 1 21 504 

11J_/ _L £i J_ UUt: 


1.18 


0.16 


1.221 


0.329 


0.03 


1.51-2.84 


HD 149143 


1.20 


0.25 


1.327 


0.053 


0.00 


1.63-3.05 


HD 75289 


1.21 


0.22 


0.466 


0.048 


0.03 


1 68-3 16 

J- . \J\J r-J* -L\J 


HD 109749 


1.21 


0.26 


0.277 


0.063 


0.00 


1.68-3.16 


HD 179949 


1.21 


0.14 


0.916 


0.044 


0.02 


1.68-3.16 


HD 86081 


1.21 


0.26 


1.49 


0.035 


0.01 


1.68-3.16 


HD 149026 


1.30 


0.36 


0.337 


0.043 


0.00 


2.31-4.34 


HD 224693 


1.33 


0.34 


0.718 


0.192 


0.04 


2.57-4.82 


Tau Boo 


1.35 


0.23 


4.126 


0.048 


0.02 


2.76-5.17 


OGLE-TR-132 


1.35 


0.43 


1.19 


0.031 


0.00 


2.76-5.17 


HD 177830 


1.46 


0.54 


1.531 


1.227 


0.10 


4.04-7.58 


Eps Ret 


1.49 


0.42 


1.556 


1.270 


0.06 


4.48-8.41 


HD 104985 


1.60 


-0.35 


6.315 


0.779 


0.03 


6.55-12.27 


1 Data from Butler ct al (2006) 


, www. exoplanet 


.eu, and references therein. 









